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Abstract 

We outline a method of deriving boost invariant hamiltonians for effective particles in quantum field 
theory The method is based on the similarity renormalization group transformation for hamiltonians in the 
canonical light-front quantization scheme. The hamiltonians are defined and calculated using creation and 
annihilation operators. The renormalization group equations are written for a sequence of unitary transfor- 
mations which gradually transform the bare canonical creation and annihilation operators of a local theory 
to the creation and annihilation operators of effective particles in an effective theory with the same dynam- 
ical content but a finite range of energy transfers due to form factors in the interaction vertices. The form 
factors result from the similarity renormalization group flow of effective hamiltonians. The regularized ini- 
tial hamiltonian and the renormalized effective hamiltonians possess seven kinematical Poincare symmetries 
specific to the light-front quantization scheme. Thus, the effective interactions can be used to describe the 
constituent dynamics in relativistically moving systems including the rest and the infinite momentum frame. 
Solutions to the general equations for the effective hamiltonians are illustrated in perturbation theory by 
second-order calculations of self-energy and two-particle interaction terms in Yukawa theory, QED and QCD. 
In Yukawa theory, one obtains the generalized Yukawa potential including its full off-encrgy-shcll extension 
and form factors in the vertices. In QED, the effective hamiltonian eigenvalue problem converges for small 
coupling constants to the Schrodinger equation but the typical relativistic ultraviolet singularities at short 
distances between constituents are regularized by the similarity form factors. In the second-order QCD 
effective hamiltonian one obtains a boost invariant logarithmically confining quark-anti-quark interaction 
term which may remain uncanceled in the non-abelian dynamics of effective quarks and gluons. 

PACS Number: ll.10.Gh 
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1. INTRODUCTION 



This paper describes a theory of boost invariant effective hamiltonians in quantum field theory. Physical 
states are assumed to be describable by solutions to the Schrodingcr equation with these hamiltonians. 

The hamiltonians are derived by integrating first-order differential equations of the similarity renormal- 
ization group. The initial condition is provided by the regularized canonical light-front hamiltonian of a 
local quantum field theory with counterterms. The hamiltonian acts in the Fock space which is constructed 
by applying the canonical bare creation operators to a vacuum state. The vacuum is annihilated by the 
corresponding annihilation operators. Due to the light-front boost symmetry, this representation of states is 
useful in a relativistic theory and the dynamics in any frame resembles dynamics in the infinite momentum 
frame. 

The similarity rcnormalization group is defined in terms of running creation and annihilation operators. 
The running operators interpolate between the bare ones in a local hamiltonian and effective ones in the 
effective hamiltonian. The effective operators are used for construction of the effective basis states in the Fock 
space. The effective hamiltonians are calculable term by term using methods of successive approximations 
and perturbation theory. One has to study the role of different terms in the effective Schrodingcr equation. 
Solutions include bound states. 

One reason for the renormalization group to play an important role in the hamiltonian approach is 
ultraviolet divergences; the initial expressions for bare hamiltonians contain divergences of local field theory 
and the divergences require renormalization. The most prominent example is the canonical hamiltonian 
of QED. Old-fashioned tree diagrams of hamiltonian perturbation theory are finite and closely reproduce 
experimental data. Renormalization problems appear when one sums over intermediate states and the 
sum diverges. The divergences correspond to diverging loop integrals in Feynman diagrams. However, the 
hamiltonian approach greatly differs from the lagrangian diagrammatic approaches. 

One apparent difference is that the sums over intermediate states involve integrals over a three-dimensional 
momentum space while the integrals in the lagrangian calculus are four-dimensional. Although a connection 
exists for finite integrals which are not sensitive to cutoffs, the connection is broken when the integrals 
diverge. In the hamiltonian approach, we introduce the effective Fock space basis and we construct hamil- 
tonians using three-dimensional regularization and renormalization procedures. In the lagrangian approach, 
one directly calculates Green functions using four-dimensional regularization and renormalization techniques. 
Equivalence of the two approaches in diverging cases remains to be shown, especially when the bound states 
are taken into account. 

Besides removing divergences, the renormalization group is useful in the hamiltonian approach because 
it introduces a hierarchy of scales. Phenomena of different scales are dealt with in a certain order. This 
enables us to solve problems involving many scales. Particle theories contain many, possibly even infinitely 
many different scales. Using the renormalization group approach, we can start from the hamiltonian of a 
basic theory that couples all degrees of freedom of all scales and we can reduce the initial hamiltonian to 
an effective renormalized hamiltonian in which the couplings between degrees of freedom of vastly different 
scales vanish. Then, the couplings between different scales are further reduced in the renormalization group 
flow to obtain the effective hamiltonian matrix which is sufficiently narrow in scale so that its spectrum of 
eigenstates can be found in practice. The similarity transformation is designed to eliminate all large changes 
of scale by the effective interactions. 

The different scales in the hamiltonian renormalization group approach are defined by different scales 
of momentum. The momentum scales are defined using the relative momenta of interacting particles. The 
definitions will be given in Section 2. Here we need to mention that the effective hamiltonian of a small 
width contains interactions which couple particles of similar energies only. The energy changes induced 
by the hamiltonian are limited by the running similarity renormalization group cutoff, denoted by A. The 
smaller is A the smaller is the energy width of the effective hamiltonian. Exact results for physical quantities 
are independent of A. 

A number of model subspaces of the Fock space need to be considered when one is solving for the spectrum 
of a field-theoretic effective hamiltonian because the full space of states is too large for computations. 
Different physical problems require different model subspaces. Working within a subspace of interest, one 
should secure that the results for physical quantities are independent of the running cutoff A. The cutoff 
independence can appear only in a certain range of cutoffs that corresponds to the physical problem and 
model subspace under consideration. However, once the cutoff independence in the finite range is achieved, 
one expects to have solved the theory in this range. 
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The following diagram will illustrate the situation. 
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H(oo,h,6, A) 
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limited 



H {X,n,S,A) 



H (\,n,6,A) 



In this diagram, the vertical arrows indicate evolution in the renormalization group parameter A which 
limits the relative energy transfers in the interaction terms. It will be explained in detail in next Sections 
how the limits are imposed. A ranges from infinity in the initial hamiltonian to a finite value in an effective 
hamiltonian. The hamiltonians depend on additional parameters n, 5 and A. 

n stands for the cutoff on the change of the particle number. It defines the limits on the numbers of 
creation and annihilation operators that can appear in a single hamiltonian term. For example, the canonical 
expressions for light-front hamiltonians in local field theories of physical interest have the number of creation 
and annihilation operators in a single term limited to 4, and the particle number cannot change by more 
than n = 2. 

S stands for the infrared cutoff. For example, it may be the lower bound on the longitudinal momentum 
carried by a particle that appears or disappears in a single interaction. 

A stands for the ultraviolet cutoff which defines the upper limit on the relative transverse momentum of 
particles which can appear or disappear in a single interaction. 

The left branch of the diagram is marked "full" because it represents the renormalization group flow 
calculated using the effective creation and annihilation operators with no restriction imposed on the space 
of states. 

The right branch of the figure is marked "limited" because it describes the renormalization group flow in 
the bare model space which is limited by parameters h, d and A. Imposing the limits is denoted by the arrow 
marked 1. For example, h can limit the number of bare particles, 5 can limit the bare particle momenta 
from below and A can limit from above free energies of the states which are taken into account. 

The initial -ff (oo, n, <5, A) contains counterterms which are constructed using the condition that physical 
results have well defined limits when the cutoffs n, S and A are relaxed. The construction of counterterms 
in perturbation theory will be discussed in detail in the next Sections. Once the counterterms remove the 
regularization dependence from the effective dynamics the arguments n, 5 and A in H(X, n, <5, A) in the 
lower left corner of the diagram are equivalent to their limiting values, n = oo, (5 = and A = oo. Thus, 
H(\,n,5, A) = H(X). One should stress that the infrared regulator <5 may still appear in the effective 
hamiltonian if there are massless particles in the theory. This is important in gauge theories. 

Accuracy of the step denoted by the arrow 1 has to be checked by relaxing the model cutoff parameters 
h, 5 and A and measuring the resulting changes in the spectrum of H(\, n, 5, A). Naturally, these cutoffs 
may have to be varied in a big range because they are introduced along the arrow 1 for bare particles. 

The effective hamiltonians at the bottom of the diagram, namely, H(X, n, S, A) and H (A, n, 5, A) with 
energy transfers limited by finite A, are connected by the arrow marked 2. This arrow denotes the procedure 
of introducing the small space cutoffs n, 5 and A which enable us to approximately solve for the spectrum 
of the effective hamiltonian H (A, n, 5, A). This time, however, the final computation cutoffs are introduced 
at the level of the effective particles, not at the level of the initial bare particles. 

The arrow 2 denotes the replacement of the whole effective hamiltonian matrix by a limited matrix. The 
procedure of obtaining the small matrix will be discussed below. The spectrum of the small matrix may be 
very close to the corresponding part of the spectrum of the full matrix because A is small (this will become 
clear later). The accuracy of the calculation must be verified by relaxing the cutoffs h, S and A and observing 
convergence of results as in the case of the arrow I and branch "limited" . But now, it is natural to expect 
that the cutoffs n, <5 and A may have to be varied only in a small range which corresponds to A. Thus, a 
finite dynamical problem to solve is defined. 

The "full" renormalization group evolution is calculable using the method described in this paper. The 
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"limited" evolution can be calculated using the matrix elements techniques introduced earlier by Glazek and 
Wilson in Refs. JlJ and || who drew on the work of Wilson ||] [Q. The matrix elements techniques were 
introduced for application to QCD. || Alternatively, one can adopt Wegner's flow equations for hamiltonian 
matrix elements in cases soluble with the energy-independent width A. || The present approach can 
be viewed as a special case of the general similarity renormalization group for hamiltonians because the 
hamiltonians we consider transform by the same unitary transformations as our creation and annihilation 
operators. However, by having defined the renormalization group transformation for the effective creation and 
annihilation operators, we remove the need to consider the model space dependence of the renormalization 
group transformation. 

The transformation we describe in this paper is partly similar to the transformation discussed by Melosh. 
|^| The important difference is that we provide a dynamical theory of the transformation in a form applicable 
to particles of different kinds. If one restricts attention to QCD, the boost invariant calculus is expected 
to help in establishing a connection between the constituent quark model, Feynman parton model, and 
perturbative quantum chromodynamics. 

Both ways in the diagram which start from the initial hamiltonian H(oo, n, 5, A) and go through the 
arrow 1 and the arrows "limited" (called branch 11) and the arrows "full" and the arrow 2 (called branch f2), 
lead to a finite H (A, h, 5, A). When calculations of some selected matrix elements of the effective hamiltonian 
are done in perturbation theory, both ways of going through the diagram are equivalent. For only a finite 
range of particle numbers and momenta can be reached in a limited number of steps of size A starting from 
the finite values selected by the external states of the matrix elements in question. 

Differences arise when one solves for the spectrum of an effective hamiltonian and when one attempts 
to vary the model space parameters h, S and A. In the "full" calculation, one obtains a single effective 
hamiltonian which one can solve in successively cnlargcablc model spaces. In the "limited" calculation, the 
model space restrictions are imposed at the beginning and they lead to an effective hamiltonian whose action 
cannot be considered in a larger model space without repeating the renormalization group calculation in the 
larger space. 

An explicit example of a difference between the two branches 11 and f2 in the diagram above is provided 
in Ref. H which discusses a Tamm-Dancoff (TD) approach analogous to the branch 11 (cf. Ref. [|l0)). In 
the TD approach, there are restrictions on the particle number which naturally lead to the sector-dependent 
counterterms as described in Ref. ||, for example, for masses. On the other hand, in the procedure of the 
branch f2 no such sector dependent counterterms arise. The present paper describes examples of sector- 
independent mass counterterms. 

Proportionality to different powers of the coupling constant helps in estimates of how important are 
different effective interaction terms and how to choose the model space. Finding the basis which can span 
a good approximation to the full solution requires trial and error studies. This general feature can be 
illustrated by the following 2x2 matrix. 



a + bg 2 gv 
gv c + dg 2 

This matrix is a model of the entire effective hamiltonian matrix calculated to second order in g including all 
couplings between all effective Fock sectors as given by the "full" calculation. Thus, we have the hamiltonian 
terms order 1, order g and order g 2 . In a perturbative calculation using matrix elements, which is focused 
on the upper sectors, one would calculate only the terms a and bg 2 . 

Assume that a and c are of the same order, and b and d are of the same order, and calculate eigenvalues of 
the model matrix neglecting terms order <? 4 and higher powers of g. For arbitrarily small g, the eigenvalues 
are given by a and terms quadratic in g. No terms linear in g arise. The quadratic corrections include 
contributions due to the term gv which couples different sectors. The role of this coupling needs to be 
estimated. The presence of dg 2 seems to be irrelevant because it couples to the upper sector through the 
off-diagonal terms order g. Hence, it seems to contribute only in order g A to the eigenvalues. 

It is well known that the above analysis is wrong in the case with degenerate diagonal matrix elements 
no matter how the degeneracy arises. For example, consider the case of a finite g such that a + bg 2 = c + dg 2 . 
The eigenvalues are equal a + bg 2 ± gv. They are linear in g instead of being quadratic, for arbitrary v. The 
lowest eigenvalue eigenstate is a superposition of the upper and the lower sector instead of being dominated 
by the upper one. In this example, the degeneracy is not visible until the term dg 2 is included in the 
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calculation. As a second example consider the case with degenerate matrix elements c and a + bg 2 and d 
not included. The simple non-degenerate perturbative expansion is again not applicable. But the addition 
of the term dg 2 can lift the degeneracy and make the simple perturbative expansion work. 

Corrections due to interaction terms such as gv may be additionally suppressed for very small A since 
the range of v in momentum variables is given by A. If A is reduced in the renormalization group flow down 
to a number on the order of some positive power of g then the resulting interaction can contribute to the 
eigenvalues in the order implied by g and A together which is higher than g 2 . In addition, the effective 
interaction v may contain small factors. For example, in the effective e + e _ -sector of positronium, the 
emission of photons is proportional to the velocity of electrons which is order a, on average. The interaction 
term gv which couples states with an additional photon, plays no role in the eigenvalue in order a 2 if the 
width A restricts energy changes to order a 2 m e i ectron (cf. 0). 

Terms such as dg 2 have been originally discussed in the light-front approach to QCD by Perry. [Q 
Heavy quarkonia are dominated by the effective QQ sector. Terms such as dg 2 in other Fock sectors may 
lift up energies of effective gluons due to the non-abelian interactions to a sufficiently high value so that the 
model hamiltonian a + bg 2 in the QQ sector alone may have eigenstates which approximate the full solution 
for heavy mesons. The important observation made by Perry in a frame dependent matrix elements 
approach using coupling coherence is that the terms bg 2 contain a logarithmically confining potential. An 
analogous boost invariant logarithmic interaction term in the Fock space in our approach will be discussed 
in Section 3. 

The above diagram and the 2x2 matrix model illustrate the structure of our similarity renormalization 
group approach to the light-front hamiltonian dynamics in quantum field theory. We summarize the steps 
here. 

The first step is the calculation of the effective hamiltonian, 

HW = Sl nAA JI(oo,M,A) S WiA . (1.1) 

S denotes the similarity transformation. Eq. (1.1) corresponds to the arrows marked full in the diagram. 
The second step is to solve the effective Schrodinger equation 

H(X) |V>) = E\1>). (1.2) 

H (A) has the same dynamical content and eigenvalues as H(oo, n, 5, A). The eigenvalue E is independent of 
the width A. 

Equation (1.2) greatly differs from the eigenvalue equation for iJ(oo, n, <5, A). The major difference is 
that the dynamics of H(X) has a limited range on the energy scale and the hamiltonian does not contain 
ultraviolet divergences. Therefore, one can attempt to solve the eigenvalue problem scale by scale. Scattering 
processes are described by the same hamiltonian. Next Sections will give examples of two fermions scattering 
in different theories. 

Solutions to Eq. (1.2) provide renormalization conditions for the finite parts of counterterms. A general 
method is necessary for reducing the full eigenvalue problem to a manageable one. This step is marked by 
the arrow 2 in the diagram. In the case of the 2x2 matrix model, this step corresponds to the calculation of 
the model space hamiltonian in the upper-left corner of the matrix. The similarity renormalization scheme 
guarantees that this step is free from ultraviolet divergences because the width A is finite. 

In the general case, one can apply the well known Bloch (l^] technique of calculating model space 
hamiltonians. Suppose we want to evaluate a model two-body hamiltonian knowing H(X) with A < to, 
where m is the effective one-body mass. We can introduce the projection operator P on the effective two- 
particle sector with a limited center-of-mass energy. We also introduce the operator R which generates the 
multi-particle and high energy components of the eigenstates from their limited mass two-body part. By 
assumption, R satisfies the conditions (1 — P)R = RP = R and PR = R(l — P) = and the equation 
(P + R — 1) H (A) (P + R) = 0. Then, the model two-body dynamics is described by the hamiltonian Q 

H 2 = {P + R^R)- 1 ' 2 (P + i? f ) H{\) (P + R) (P + R ] R)- 1/2 . (1.3) 

The same approach can be used for larger model spaces. The model space is characterized by the parameters 
h, S and A in the diagram. So, the operation R depends on these parameters. But the resulting spectrum 
in the range of interest should not depend on the model space boundary when the width A is small and the 
model space contains the dynamically dominant basis states in the selected range of scales. The heuristic Eq. 
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(1.3) can be applied in perturbation theory in the effective interaction even for sizable coupling constants 
since the effective interaction strength is considerably reduced by the similarity factors. 

The scheme outlined above is still prone to the infrared rcgularization dependence for massless particles. 
This is particularly important in gauge theories. However, the effective hamiltonian dynamics is expected 
to lead to infrared convergent results for gauge invariant quantities. There is also a possibility that new 
effective interactions are generated from the infrared region and they bring in effects normally associated with 
a nontrivial vacuum state. || [fl4[ |l5[| We shall make comments on the issue of long distance phenomena in 
the present approach in Section 2.b where we describe the range of scales involved in the theory. The reader 
should refer to j|3], [p^[ , Q and jl(| for discussions of the ground state, spontaneous symmetry breaking 
and zero-modes problems. 

The remaining part of this paper is organized as follows. Section 2 describes the hamiltonian formalism 
in three subsections. Namely, Section 2. a introduces the similarity renormalization group equations and 
describes methods of solution, Section 2.b describes regularization factors, 2.c deals with renormalization 
conditions. Section 3 contains examples of lowest order applications of the formalism. Our derivation of 
the generalized Yukawa potential is given in Section 3. a. Section 3.b describes the Schrodinger equation 
for positronium in QED. Section 3.c discusses a confining term for constituent quarks in QCD. Section 4 
concludes the paper. The list of references is focused on the similarity renormalization group approach to 
hamiltonian dynamics in the light-front Fock space. The reader should be aware of this limitation. Examples 
of other approaches to quantum field theory in the light-front form of dynamics can be found in Ref . |lG| . 

2. EFFECTIVE HAMILTONIANS 

This Section is divided into three parts. The first part describes our method of calculating effective 
hamiltonians in the Fock space. The second part presents our regularization scheme for initial hamiltonians. 
The last part discusses renormalization conditions and the effective eigenvalue problem. 

2. a Similarity transformation 

We construct a family of effective hamiltonians in the light-front Fock space. The family is parameterized 
by a scale parameter A which ranges from infinity to a finite value. A limits energy transfers in the interaction 
terms. 

The hamiltonians are built of sums of ordered products of creation and annihilation operators. The 
hamiltonian labeled by A is expressed in terms of creation and annihilation operators which correspond to 
A. We commonly denote these operators by q x . In addition, the creation and annihilation operators carry 
labels of quantum numbers such as momentum, spin, flavor or color. We will not indicate those numbers in 
the initial presentation, unless it is necessary. 

All hamiltonians in the family are assumed to be equal. Thus, 

H Xl (q Xl ) = H X2 (q X2 ). (2.1) 

For A = oo, the hamiltonians Hoc are expressed in terms of operators creating and annihilating bare 
particles, qoo. Hamiltonians can be constructed from the canonical field theoretic expressions for the 
energy-momentum density tensors. 

Unfortunately, expressions for H^, in local field theories are divergent. They need to be regularized by 
introducing a bare ultraviolet cutoff which we shall denote by e. The ultraviolet cutoff A from the previous 
Section corresponds to A 2 /e where A is an arbitrary finite constant which carries the necessary dimension 
of a mass. The limit of removing the bare ultraviolet cutoff will correspond to e — > 0. 

H x=oa — H e for all values of e. For the limit e — > to exist the hamiltonians must include a number 
of additional terms (called counterterms) whose structure will be determined later. 

Haa may include an infrared regulator, generically denoted by S. For example, this is required in QED 
with massless photons and in QCD with massless gluons. S — > when the infrared regularization is removed. 
The parameter 8 is indicated explicitly if needed. 

Our key assumption is that the particle degrees of freedom for all different scales A are unitarily equivalent 
to the bare particle degrees of freedom: 



G 



(2.2) 



This assumption says that the quantum numbers of bare and effective particles are the same for all values of A. 
The following examples explain the origin of this assumption. (1) Constituent quarks have the same quantum 
numbers as current quarks. (2) We use the same quantum numbers for photons and electrons independently 
of the kind of processes we consider in QED or in related effective theories such as the nonrelativistic 
Schrodinger equation with Coulomb potentials between charges. (3) Pions and nucleons in nuclear physics 
have the same quantum numbers quite independently of what kind of interactions, pion-nucleon vertex form 
factors or other dynamical assumptions one uses. 

It follows from Eq. (2.2) that creation and annihilation operators for Ai and A2 are unitarily equivalent 
and connected by transformations of the form Ux^U^ . The transformations U\ 1 or U\ 2 will depend on the 

bare cutoffs but the transformation U\Jj\ for finite Ai and A2 will have to be finite in the limit e — > 0. 

The transformation U\ is defined indirectly through a differential equation of the type used in the 
similarity renormalization scheme for hamiltonians of Glazek and Wilson |IJ Q . That scheme was originally 
developed for application to QCD ||. Hamiltonians with labels Ai and A2 are connected by integration of 
the differential equation from Ai to A2 . Our guiding principle in writing the differential equation for effective 
hamiltonians is that the resulting interactions between effective particles with considerably different scales of 
relative momenta are suppressed. 

This principle has its origin in the following examples. (1) Emission and absorption of short wavelength 
photons are not essential in the formation of atoms. (2) Emission and absorption of hard pions by nucleons is 
not important in nonrelativistic nuclear physics. (3) Constituent quarks have moderate momenta and their 
effective dynamics seem to be independent of the very hard gluon emissions. (4) High momentum transfer 
phenomena are independent of the small momentum transfer effects such as binding. A standard way of 
achieving this kind of picture in theoretical models is to include form factors in the interaction vertices. The 
form factors quickly tend to zero when momenta change by more then the size of a specific cutoff parameter. 

The cutoff parameter in the form factors sets the scale for allowed changes of momenta. It determines 
the range or width of the interaction in momentum space. That width is the origin of our scale A which 
labels renormalized effective hamiltonians. Our similarity factors are analogous to the vertex form factors 
which are commonly used in nonlocal models (see also Ref. [0). The large momentum transfer dynamics 
is integrated out through the similarity renormalization group equation. 

Boost invariance requires that the individual momenta of effective particles are not restricted because 
boosts change those unlimitedly. The hamiltonian width restricts only relative momenta of effective particles. 
Also, the larger is a relative momentum the larger change is generated by a boost. Therefore, when the free 
energy of interacting particles in their center-of-mass frame (i.e. the free light-front invariant mass) is much 
larger than A the immediate change of energy due to interaction is limited by the large energy itself instead 
of A. At the same time, this condition takes care of the property of wave mechanics that strong interference 
occurs between waves of similar wavelengths within a range of wavelengths on the order of the wavelengths 
themselves. || 

In our construction, strong dynamical interference effects for states of similar free energies are made 
not to contribute in the derivation of effective hamiltonians. For example, the similarity transformation is 
constructed in such a way that only large energy denominators can appear in the perturbative calculations 
of effective hamiltonians and small denominators are excluded. Namely, only large free energy changes 

are integrated out. In this approach, calculations of strong coherence effects for nearly degenerate states 
are relegated to a later step of solving for eigenstates of the effective hamiltonian. That step may be non- 
perturbative. For example, the Coulomb potential of QED is formally of the first order in a and leads to a 
variety of bound atomic structures beyond perturbation theory. 

Our differential equations require a separation of the changes in creation and annihilation operators from 
changes in coefficients in front of products of the operators. In order to define this separation we assume 
that terms with a large number of the operators in a product do not dominate or mediate the effective 
dynamics of interest. If the latter assumption turns out to be invalid our formalism merely provides a way 
to approach the resulting problems. The comment due here is that if the dynamics leads to spontaneous 
symmetry breaking, or condensates, we will have a well defined renormalized hamiltonian theory to study 
those phenomena in the desired detail, cf. Refs. Q, |^5| and ||. 

The unitary equivalence of creation and annihilation operators at the scale A and at the infinite scale, 
i.e. those appearing in = H e , together with the equality of hamiltonians at all scales imply that 
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H X (q X ) = U x H x (q 0Q )U{ = H 0O {q 0O ). (2.3) 
We denote H x (q OCl ) = TL X and obtain 

H x = UlH^Ux. (2.4) 

Thus, the effective hamiltonian H\ is obtained from the hamiltonian 7i x by replacing creation and anni- 
hilation operators for bare particles by creation and annihilation operators for effective particles with the 
same quantum numbers. The bare creation and annihilation operators are independent of A. One calculates 
A-dependent coefficients in front of the products of q^ in TL X . 
The differential equation for TL X is Q 

^■H X = [H X ,T X ], (2.5) 



where 



7i\ has the following structure 



T> = UI±U X . (2.6) 



Hx = F x [Q x ]. (2.7) 

■F\[£?a] denotes the similarity form factors in 7i x to be described below. Using the unitary equivalence, we 
also have 

H x (q x ) =F x [G x (q x )}, (2.8) 

where 

G x (q x ) = U x G x ul (2.9) 

A similar relation holds for T x (q x ) and T x since the latter is expressed in terms of q^. 

The operation F x acts on the operator Q x equal to a superposition of terms each of which is an ordered 
product of creation and annihilation operators. The ordering is arbitrary but needs to be determined. We 
adopt the order from left to right of creators of fermions, creators of bosons, creators of anti-fermions, 
annihilators of anti-fermions, annihilators of bosons, annihilators of fermions. At least two operators must 
appear in a product and at least one creation and one annihilation operator must appear. No product 
contains only creation or only annihilation operators. This is a special property of light-front hamiltonians. 
Hamiltonians in other forms of dynamics do not have this property and lead to the necessity of solving 
the ground state formation problem before other states can be considered because the pure creation or 
annihilation terms produce disconnected vacuum dynamics. 

The operator Q x is divided into two parts, Gi X and G2\- Gi\ is a superposition of all terms of the form 
a' a for a equal q^ of any kind. In principle, one could also include in Gi X some chosen terms with a larger 
number of creation and annihilation operators, e.g. terms containing two creation and two annihilation 
operators. However, plane-wave Fock space basis states are not eigenstates of relevant operators of such 
type and we limit Gix to terms a^a to avoid the difficulty in present calculations. 

Gi X becomes the effective free part of G x , denoted G\ X , after q^ is replaced by q x . The effective free 
hamiltonian part Hi X is equal to Gia because G\ X is not changed by the operation F. Eigenvalues of G\ X 
are called free energies. 

The remaining part G2\ = Gx — Gix gives the interaction part of the effective hamiltonian H x . One 
replaces goo by q x and obtains G*2A ■ Then, one applies the operation F x which inserts the vertex form factors 
defined as follows. 

Let the momentum labels of all creation operators in a single product in an interaction term be k% , k2, kj 
and the momentum labels of all annihilation operators be k[, k' 2 , kj. Each momentum has three compo- 
nents, k + ranging from to oo and two transverse components k 1 - = (fc 1 ,/; 2 ), both ranging from — oo to 
+oo. The z-axis is distinguished by our choice of the light-front. The sum of momentum labels of creation 
operators, 2i=i^ii equals the sum of momentum labels of annihilation operators, Ylj=i^j- We denote 
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these sums by P = (P + ,P ± ), (P + is positive). Each k + or k' + is a positive fraction of P + ; Xi = k^/P + , 
1 > Xi > and x'j = k'+ /P+, 1 > x'j > 0. We have YlLi x i = E/=i 4 = l - We also denne 

4 = fc^-a;„P x (2.10) 

for all momenta in the hamiltonian term. Ei=i K i~ = Ej=i = 0- 

These momentum variables appear standard but the way they are used here is not. Namely, P is usually 
not equal to a total momentum of a physical state. It characterizes the interaction term whose action 
redistributes P from the set of momenta of the annihilated particles to the set of momenta of the created 
particles. 

Thus, each term in the hamiltonian is characterized by P and two sets of variables, Xj = {(x i: 
for creation operators and X'j = {Ocj; ft'j )K'=i for annihilation operators. For example, in a product of 
two creation operators and one annihilation operator we have x\ — x, X2 = 1 — x and x[ — 1. Also, 

= = k and re'j =0. P can be arbitrary and the term in question replaces one particle of 

momentum P by two particles of momenta xP + k and (1 — x)P — k for + and _L components, k + = 0. It 
is convenient to speak of P as a parent momentum and about the individual particle momenta as daughter 
momenta. The parent momentum in a hamiltonian term equals one half of the sum of momenta labeling 
all creation and annihilation operators in the term. Each daughter particle carries a fraction of the parent 
momentum. The parent momentum may be carried by one or more particles. 

The operation F\ acting on a product of creation and annihilation operators produces 



n °l n a ^ 

i=i j=i 



MX^X'jUalUa*,. (2.11) 

i=i j=i 



The function f\(Xi, X'j) is a suitable function which represents our physical intuition about form factors. 
The arguments of fx are invariant with respect to seven kinematical Poincare transformations of the light- 
front frame. This feature results in the boost symmetry of our theory. We impose three conditions on the 
function f\. 

The first condition is that fx is expressible through the eigenvalues of G\\ corresponding to the sets 
Xi and X'j so that fx equals 1 for small differences between the eigenvalues and quickly goes to zero 
when the differences become large. This is the basic condition of the similarity renormalization scheme for 
hamiltonians.0 [0] The width of fx is set by A. One can consider functions fx which depend on Xj and 
X'j in a more general way than through the eigenvalues of the free hamiltonian but that option will not be 
investigated here. 

The first condition defines the effective nature of the hamiltonian labeled by A. Namely, the effective 
particle states which are separated by the free energy gap which is much larger than A arc not directly coupled 
by the interactions. In other words, A limits the free energy changes induced by the effective interaction. 
Moreover, as a consequence of fx ~ 1 for similar energies, 1 — fx is close to zero for the similar energies and 
it vanishes proportionally to a power of the energy difference. The higher is the power the smaller is the 
role of states of similar energies in the calculation of the effective hamiltonian. This will be explained later. 
Consequently, the higher is the power the smaller is the role of non-perturbative phenomena due to energy 
changes below the scale A in the calculation of the effective hamiltonian. 

Thus, there is a chance for the full hamiltonian diagonalization process to be divided into two parts: a 
perturbative calculation of the effective renormalized hamiltonian and a non-perturbative diagonalization of 
that effective hamiltonian. This is our factorization hypothesis in the hamiltonian approach. 

The second condition is that both, 1 — fx and dfx/dX, must vanish faster than linearly in the free energy 
difference. This condition is required to exclude the small energy denominators in perturbation theory and 
will be explained below. The second condition implies that 1 — fx vanishes as at least second power of the 
energy difference near zero. 

The third condition is defined by saying that multi-particle interactions (especially interactions that 
change the number of effective particles by many) should not be important in the effective hamiltonian 
dynamics which is characterized by changes of energies below the scale A. This may be possible if fx as 
a function of the daughter variables approximates the shape of one particle irreducible vertices which is 
characteristic to the theory under consideration. Structure of Gx depends on the choice of the function fx. ■ 
Some choices will lead to more complicated effective interactions than others. The best choices for the most 
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efficient description of physical phenomena at some scale A are such that the effective particles interact in 
a way that is most easy to understand and which can be parametrized with the least possible number of 
parameters over the range of scales of physical interest. One can conceive variational estimates for the best 
choices of fx that minimize complexity of the effective hamiltonians. For example, it is clear that creation of 
effective particles will be suppressed when the width A becomes comparable to the effective masses of those 
particles. 

To satisfy the first condition above in a boost invariant way we define a boost invariant gap between free 
energy eigenvalues for effective particles which is to be compared with the running cutoff parameter A. The 
light-front quantization scheme does not explicitly preserve rotational symmetry. Nevertheless, it has been 
shown that if counterterms provide enough freedom through their finite parts and multi-particle effects are 



suppressed one can obtain rotationally invariant results. JL8I |19 |l 



X'j are 



The free energy eigenvalues relevant to a particular hamiltonian term with daughter variables Xj and 

+ mf(A) _ P ±2 +Mj 



E 



and 



(2.12) 



(2.13) 



where 



„X2 



(=1 



and 



/2 J K >f 2 + m 2 J\) , 
M' 2 j = J2 x , ■ ( 2 - 15 ) 



The individual effective particle masses are allowed to depend on the effective hamiltonian width parameter 
A. We define the mass difference for a hamiltonian term to be 

AM 2 =M' 2 j-Mj, (2.16) 

and the mass sum to be 

EM 2 = M' 2 j + Mj. (2.17) 

To be specific, we define details of the function /a(Xj,Xj) introducing a parameter z\. Following the 
similarity renormalization scheme [[l] [0, z\ can be chosen in such a way that z\ is close to zero for AA4 2 
small in comparison to A 2 or T,A4 2 and \z\\ is close to 1 for AA4 2 large in comparison to A 2 or comparable 
to T,A4 2 . For example, 

AM 2 

* = ^Ta^' (2 - 18) 

The definition includes HM 2 ) to ease estimates in high order perturbation theory, especially in the analysis 
of overlapping divergences. [Q The new feature here is that the introduction of TjM 2 does not violate the 
light-front boost invariance and basic cluster decomposition properties. f\(Xi, X'j) is defined for the purpose 
of this article to be a function of z 2 , n > 1, which is analytic in the vicinity of the interval [0, 1] on the real 
axis, equals 1 for z\ — and quickly approaches for z\ ~ 1; 

f x (X I ,X' J )=f(zT). (2.19) 

For example, 
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/(«) 



u(l - Up) 
u (l - u) 



(2.20) 



where 1 > it > and m > 1. The larger the exponent m the closer f(u) approaches 6(uq — u) for < u < 1. 
Eq. (2.20) concludes our definition of the operation F\. 

The smallest possible value of TjM 2 in Eq. (2.18) is ^2 i=1 m iW + 2j=i m jW ■ Thus, is small 

for small positive A 2 when AM 2 is small in comparison to particle masses. Therefore, ito must be much 
smaller than 1 to force AM 2 to be small in comparison to XLM 2 when A 2 is small. One can also force AM 2 
to be small in comparison to the particle masses by making A 2 negative so that it subtracts from Y>M 2 its 
minimal value. Then, the mass difference is compared to the sum of kinetic energies due to the relative 
motion only. It is also useful to limit the small mass differences by choosing an infinitesimally small uo and 
introducing A 2 = m 1/>2 A 2 . Then, |AA1 2 | < A 2 in the ^-function limit. In this case, the band-diagonal 
hamiltonian width becomes independent of the mass sum for as long as the latter is small in comparison to 
A 2 . 

The infinitesimal transformation T\ in Eq. (2.5) is defined as follows. Eq. (2.5) is rewritten using Eqs. 
(2.8) and (2.9), into the form 



W = f'Q + SQ' = T] + [fG 2 ,T] 



(2.21) 



The prime denotes differentiation with respect to A. We have simplified the notation of F\{G\] to fQ. Three 
universal relations fQi — Q\ and (1 — f)Q\ = f'Qi = are then used without saying. 

Equation (2.21) involves two unknowns, Q and T. Additional condition is required to define T. One 
recalls that if the interaction is absent, i.e. when Q2 — 0, then no evolution with A may appear. Therefore, 
in the limit of negligible interactions, both Q' — and T = . Q' should differ from zero if and only if 
the interactions are important. The first term on the right-hand side is order T since Q\ contains terms 
independent of interactions. The second term on the right-hand side is at least of second order in interactions. 
The first term can be used as a seed for defining T through a series of powers of the interaction. 

We associate the derivative of Q with the second term on the right-hand side. The first term on the 
right-hand side and a part of the second term which is left after the derivative of Q is defined, together 
determine T. T is defined through the commutator [Gi,T] using a curly bracket notation. We write 



when 



A = {B} Ql 



[A, Gx] = B. 



(2.22) 



(2.23) 



Subscripts of such curly brackets are often omitted in later discussion. Suppose B contains a term which 
involves a product 

7 ./ 

1=1 j=i 

Then, {B}g lx contains the same product (as a part of the same expression) with an additional factor equal 



£ k'f 2 +rn 2 (\) ^ fc ±2 + m 2 (A) 



kt 



(2.25) 



The sums of individual energies satisfy Eqs. (2.12) and (2.13), respectively, and the factor (2.25) equals 



AM' 



(2.26) 



where P + is the parent momentum for the product under consideration and the mass difference is defined 
in Eq. (2.16). All terms in the operator B are multiplied by the corresponding factors. 
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The factor (2.25) explodes to infinity when the denominator approaches zero. Hence, for the operator 
A to be well defined, the coefficients of products of the form (2.24) in the operator B must vanish at least 
as fast as the energy denominator itself when the denominator approaches zero. Therefore, our definition of 
[Gi,T] is given in terms of an operator which has such property. Eq. (2.21) is split into two equations as 
follows. 



!9' = f[f92,T] 



(2.27) 



[T,gi] = (i-f)[fg 2 ,T]-fg. 



(2.28) 



Our second condition introduced below Eq. (2.11) on the functions 1— f\{Xi, X'j) and f' x (Xj, X'j) guarantees 
that T is well defined and tends to zero in the region of vanishing energy denominators because the right-hand 
side of Eq. (2.28) vanishes at least as fast as the first power of the energy differences. 

Equation (2.27) is a first order differential equation. One has to provide an initial condition to define a 
theory The initial conditions are set in this paper by canonical light-front hamiltonians plus counterterms. 
The latter are determined from the condition that the effective hamiltonians have well defined limits when 
the bare cutoff is removed. In other words, one has to find the class of initial conditions at A = oo which 
imply e-independent H\s for all finite A's when e — > 0. 

A general iterative procedure for calculating the effective hamiltonians is analogous to the one from 
Refs. [Q and j|. However, instead of iterating two coupled equations for H\ and T\ we base iteration on 
an equivalent single equation for H\ with an explicit solution for T\ already built in. Simple algebra and 
substitution of Eq. (2.28) into Eq. (2.27), lead to 



d a 
~d~X 



fxG: 



2A. 



d , 



h)Q: 



2 A 



(2.29) 



Equation (2.29) drives the renormalization group formalism in this paper. Note that the right-hand side is 
given in terms of a commutator. Therefore, the effective renormalized hamiltonians contain only connected 
interactions. This is essential for cluster decomposition properties of the effective hamiltonians. |20[ 
Equation (2.29) is of the form 



dX 



Qx = T x [g x ]. 



(2.30) 



The right hand side contains terms which are bilinear in the effective interaction strength. The initial 
condition for Eq. (2.29), or (2.30), is given at A = oo: Gx=oo — Qe- 

With accuracy to the first order in powers of the interaction strength, Q\ is independent of A and Q t 
is equal to the initial regularized hamiltonian expression one intends to study, denoted by . In this 



initial approximation, 7i, 



(0) _ f(0)g(0) 



where Q x °^ — and f^ is the similarity factor calculated using 



eigenvalues of G$ . Ti.^ forms our first approximation to the similarity renormalization group trajectory of 
operators Q\ parametrized by A. 

Eq. (2.30) can then be written in the iterative form for successive approximations to the trajectory Q\ . 
Namely, 



This is an abbreviated notation for 



(n+l) 



(2.31) 



d r (n+l) 

dX 



/J? denotes a function of z x n> expressed through eigenvalues of Gi\, such as in Eqs. (2.19) and (2.20). The 
initial condition is set by Q, 



» 



f {n) r (n) \ a f(n)\ n (n) 

h y 2 x »i dA u /a )y 2 x 

(n) 



(2.32) 



(n+l) 



Qz +1 \ Thus, the solution is 



?("+!) 



rp(n) \g(n) 



(2.33) 
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Goo contains the e-regulated canonical hamiltonian terms and counterterms. The counterterms remove 
the part of the integral in Eq. (2.33) which diverges for finite A when e — ► 0. Matrix elements of the 
hamiltonian of the effective theory are required to have a limit when e is made very small. The condition 
that the necessary Goo exists is the hamiltonian version of renormalizability. It does not require the number 
of counterterms to be finite, although a finite number has the clear advantage of simplicity. 

The part of the integrand in Eq. (2.33) which leads to the divergence is denoted by m n) [Gs n) ] 



and the remaining part by Ts n \Gs^] ■ Ge n+1> contains the initial regulated hamiltonian terms and 

J conv 

counterterms. The counterterms in Ge n+1 ^ are discovered from inspection of F^ n+1 ^ [G^™ +1 ' ) ] dependence on 
e when e — > in the absence of counterterms. 

Note that Fj n+1) [G^ +1) ] = G[™ +1) and it is not necessary to know F^ n+1) to calculate G^ +1) . One 

calculates F^ n+1 ^ after G^ +1 ^ is made independent of e when e — > 0. 

The diverging dependence on e when e — > 0, is typically of the form e _1 or log e times operator coefficients. 
The operator coefficients can be found by integrating the diverging part of the integrand from some arbitrary 
finite value of A, say A , to infinity. The divergence originates from the upper limit of the integration and it 
is independent of A . The remaining finite part of the integral is sensitive to the lower limit of integration 
and depends on A . The counterterm does not depend on Ao but it contains an arbitrary finite part which 
emerges in the following way. 

The counterterm subtracts the diverging part of the integral. But subtracting terms with diverging 
functions of e times known operators does not tell us what finite parts times the same operators to keep. 
Thus, one needs to add arbitrary finite parts to the numbers 1/e and loge in the counterterms. These finite 
parts are unknown theoretically and have to be fitted to data. In particular, observed symmetries may 
impose powerful constraints on the finite parts. 

The diverging part of the integrand is such that the lower limit of its integration produces the same 
operator structure as the upper limit but the diverging numbers such as e _1 or log e from the upper limit are 
replaced by finite numbers at the lower limit. Those finite numbers depend on Ao but, once they arc replaced 
by the required unknown finite parts, one obtains a valid expression for the counterterm. The replacement 
is achieved by adding to the integral the same operators multiplied by the numbers which are equal to the 
differences between the unknown numbers and the numbers resulting from the lower limit of the integration. 
Thus, the unknown numbers we need to add to the integral of the diverging part of the integrand from Ao 
to infinity depend on Ao. One can write the Ao-independent Q^ n+1 ^ as G^unite^®) Ja° ^"'[^"'l 

L J div 

The free finite parts of the counterterms are contained in G^jinitei^o) an d one can fit them to data using 
predictions obtained from effective hamiltonians at some convenient scales A. 

More than one scale A may become necessary for accurate determination of the free parameters when 
their values have to be of considerably different orders of magnitude and require knowledge of physical 
phenomena at different scales. In the present work a single scale A = Ao is sufficient for practical calculations. 
The renormalization conditions are set using H\ . One may also consider renormalization conditions for 
parameters in H\ g which are set using another effective hamiltonian at some nearby scale Ai ^ Xq- This will 
be illustrated in the next Section. 

The complete recursion including construction of counterterms in Eq. (2.33) is given by 



div 









T (n) [g ( n)] 




[ ds 


" T W [£(»)]" 


ds 




h 




div Jx 




conv 



si n+1) =0e ( XL(Ao)+ / '/* Trier: , -/ «/* ir?H n.-.u) 

In the limit n — > oo, if the limit exists, one obtains 

n X rOC 

Gx=Gefinite(Xo)+ ds[T s [G s }] dlv - ds [T s [G s ]] conv . (2.35) 

H\ is obtained from Eq. (2.35) through the replacement of by q\ (to obtain G\) and action of F\ on 
G A . 

Perturbative calculations of renormalized effective hamiltonians are based on the observation that the 
rate of change of G\ with A can be expanded in a power series in the effective interaction G2X at the same 
running scale A. This is obtained by repeated application of Eq. (2.29). One rewrites Eq. (2.29) as 
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dX 



[f\G2\, {- f'\G2\} g lx ] + /a!?2A, \ (1 - h)~7T 




(2.36) 



Then, one replaces G'\ in the last term on the right-hand side of Eq. (2.36) by the preceding terms. Two 
successive substitutions produce an expression for Q' x with four explicit powers of the effective interactions 
and the remaining terms are of higher order [note that (1 — f\)G[\ = 0]. 



We have omitted subscripts 2, A and Q\\ on the right-hand side. All the subscripts appear in the same 
pattern as in Eq. (2.36). Correspondingly, the infinite chain of substitutions produces an expression ordered 
by explicit powers of the effective interactions, to infinity. 



The round bracket raised to the n-th power means n consecutive repetitions of the symbols from within the 
round bracket. The subscripts are omitted for clarity as in Eq. (2.37). 

The above expansion in powers of the effective interactions provides a systematic order by order algorithm 
for building an expression for the effective hamiltonian. The energy denominators and functions fx arc 
calculated using eigenvalues of Gix- Therefore, in Eqs. (2.37) and (2.38), the derivatives of the function 
fx contain two kinds of contributions: those resulting from differentiating the explicit A dependence in the 
arguments z\ (for example, A 2 in Eq. (2.18) ), and those resulting from differentiating the free energy 
eigenvalues (for example, A.M 2 in Eq. (2.18)). Since the free energy terms include effective masses which 
depend on the width A, the derivatives of the effective masses appear in the equations on the right-hand 
side. Moving them to the left-hand side leads to coupled nonlinear differential equations for the effective 
hamiltonians. 

The general iterative approach in Eq. (2.34) or the expansion in Eq. (2.38), can be analysed using 
expansion in the running coupling constants. One can divide Q\\ into two parts: one which is independent 
of the coupling constants and another one which vanishes when the coupling constants are put equal to zero. 
The parts depending on the coupling constants are moved over to Q2X and treated as an interaction. After 
Qix is reduced to the part which is independent of the interactions, the derivatives of fx in Eqs. (2.36) 
to (2.38) do not introduce additional powers of the interaction strength and the series is strictly ordered 
in powers of the interactions according to their explicit appearance in the formula (2.38). This series then 
provides the perturbative expansion in terms of the running coupling constants. 

The simplest case of the perturbative expansion involves a single coupling constant at a single scale. 
Firstly, one expands the renormalization group equations into a series of terms ordered by powers of the 
bare coupling g . Secondly, one evaluates the effective coupling g\ at the chosen scale Ai as a power series 
in the bare coupling. Thirdly, the latter series is inverted and the bare coupling is expressed as a series in 
the effective coupling g\. Then, one can pursue perturbative calculations in terms of the effective coupling. 
In particular, one can reduce the hamiltonian width to A2 < Ai and calculate g2 as a series in g\. Such steps 
can be repeated. For example, one can reduce the width in each step by a factor 2. || j| N steps will reduce 
the width by the factor 2~ N . This way one can build the renormalization group flow indicated by the chains 
of small arrows in the diagram discussed in Section 1. If many coupling constants appear but they can be 
reduced to functions of a finite set of independent running coupling constants the finite set determines the 
theory and one speaks of coupling coherence. |^l| 

2.b Regularization 

A canonical bare hamiltonian obtained from a local field theory is divergent. This Section describes 
how the ultraviolet singularities in the canonical hamiltonian are regularized with the bare cutoff e. We 
also introduce infrared regularization. Our presentation is ordered as follows. First, we briefly explain 
connection between the ultraviolet and infrared regularizations in light-front dynamics. Then, we proceed 



^ = [fg, {-fg}] + [fg, {(1 - /) [fg, {-fg}}}} 

+ [fg, {(1 - /) [fg, {a - /) [fg, {-mm}] + °(s 5 ). 



(2.37) 




(2.38) 
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with definitions of the canonical hamiltonian terms. For that purpose, we have to discuss the fundamental 
set of scales in the hamiltonian approach and explain the role of lagrangian densities for classical fields in 
the construction. Then, we describe details of the ultraviolet and, subsequently, infrared regularizations. 

The ultraviolet and infrared regularizations are connected through masses. The infrared structure is 
influenced by the masses in the initial hamiltonian ff £ (0) . iJ e (0) carries the superscript to indicate that it is 
the initial hamiltonian which does not yet include counterterms. If the ultraviolet counterterms change the 
masses the infrared behavior is changed too. 

An initial mass value is generically denoted by m{°\ A light-front energy, p^, of a free particle with a 
four-momentum p m — {p + ,p^ : Pm = (p ±2 + ™ 2 )/p + ), tends to infinity when m 2 > and p + tends to zero. 
But p~ may be finite or even approach zero in this limit if m 2 = and p approaches zero too. The limit of 
small momentum p + is always a high-energy limit when m 2 > 0. But it ceases to be the high-energy limit for 
the small transverse momenta if m 2 — > 0. Thus, m e ^ in the initial hamiltonian is capable of switching from 
the high-energy regime in the longitudinal direction to the low-energy one when we take the limit m e ^ — > 0. 
Conversely, introducing masses turns the infrared low-energy regime into the high-energy regime. 

We begin the construction of with an enumeration of momentum scales. We distinguish scales 
related to the boundary conditions for fields at spatial infinity, small momentum cutoffs, phenomenological 
parameters and large momentum cutoffs. 

The bare hamiltonian is defined in terms of the operators goo. The bare quantum fields are built 



from these operators with plane-wave coefficients. 20 The initial basis in the Fock space is built from the 
vacuum state |0) using q^. Fermion, anti-fermion and boson creation and annihilation operators are denoted 
by &t , d\ , b, d and a, respectively. For example, \ka >= frj^lO) denotes a state of one bare fermion of 
momentum k = (fc + , k^). Spin z-axis projection, flavor, color or other quantum numbers, are denoted by a 
common symbol a. The momentum variable in the subscript is distinguished in order to describe the scales 
involved in the definition of the hamiltonian. The order of scales in momentum space is reverse to the order 
of scales in position space. 

The largest scale in the position space is the quantization volume. In other words, the momenta can 
be thought of as discrete when necessary. But we insure by our choices of scales that the granulation in 
momentum space is never noticeable and the quantization volume is effectively infinite for all our purposes. 
Thus, we universally adopt continuous notation for momentum variables. 

Potentially related to the boundary conditions, hypotheses about zero-modes and spontaneously broken 
symmetry in light-front quantum field theories were recently discussed in Rcf. Q which quotes important 
earlier literature on the subject. Basically, one may expect that new terms emerge in the effective hamil- 
tonians and the new terms account for the large scale dynamical effects. Susskind and co-workers have 
proposed a way to think about the wee parton dynamics in a model. |L4| Ref. |l5| describes QCD sum rules 
technique using the notion of vacuum condensates in the light-front scheme. However, the original quantum 
dynamics of the vacuum formation and spontaneous symmetry breaking are not yet understood and cannot 
be discussed further here. 

The next smaller size in position space is the inverse of the infrared regulator. Two types of the infrared 
regulator appear. The first is a lower bound, denoted by <5, on the parent -(—momentum fraction that can 
be carried by an operator in an interaction term. The second is a mass parameter fj,g. fj,$ appears as the 
mass parameter mf in the initial hamiltonian He°K It is introduced for massless bare particles, fig cuts 
off the small longitudinal momentum region at a small scale order /x 2 /A 2 , where A 2 is the invariant mass 
upper bound. When the infrared regularization is being removed, S or are sent to zero but their inverse 
is always kept negligible in comparison with the quantization volume scale. 

The next smaller scale in position space is set by the size of the volume used for preparation of incoming 
and detection of outgoing particles (including bound states) and the corresponding time scale. Physics is 
contained within this scale and observables are allowed to depend on this scale since the preparation and 
detection of states is a part of a physical process. 

The order of magnitude of momenta larger than the experimental wave packet widths are characterized 
in terms of three different scales, (1) masses of particles, (2) the width of the effective hamiltonian (i.e. A), 
and (3) the bare cutoff scale e^ 1 . When solving for the hamiltonian spectrum, a new scale may emerge 
dynamically, determined by the effective coupling constants, masses and width of the effective hamiltonian. 
Scale invariance at large momenta may be violated through a dimensional transmutation even if all mass 
scales are negligible in comparison to the momenta and A. 
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The width A ranges from oo to convenient finite values. Description of physical phenomena involving 
energy-momentum transfers of the order of k require A to be larger than k. It is also useful to use A not 
too large in comparison to k in order to avoid too much detail in the dynamics. For example, useful values 
of A in nonrelativistic systems are smaller than effective masses. In QED, the convenient A in the hydrogen 
description is much larger than the binding energy and much smaller than the electron mass. 

The bare cutoff scale e _1 is the largest momentum scale in the theory. The formal limit A — > oo is used 
only to remove A dependence from the hamiltonian regulated by e. In other words, no A dependence appears 
in the hamiltonians with A larger than the scale implied by e" 1 . No physical quantity depends on e when 
e — * 0. The similarity renormalization scheme for hamiltonians is built to achieve this goal to all orders in 
perturbation theory (cf. |jj Q). 

We proceed to the explicit construction of simplest terms in the hamiltonian . Details of counterterms 
are not known from the outset. Light-front power counting rules are helpful in determining possible 
structures of the counterterms but more is required in practice. The similarity renormalization group provides 
the required details. 

All starting hamiltonians in quantum field theories contain a free part which we denote by G[ . The free 
part for fermions and bosons has the form 



k 



k^+^ 2 f 

rr a ka a k « 



(2.39) 



We adopt the following conventions. Summation over a denotes a sum over all quantum numbers except the 
momentum. 



(2.40) 



The creation and annihilation operators in Eq. (2.39) are the bare ones denoted in Section 2. a by goo. They 
satisfy standard commutation or anti-commutation relations 



aka,a k , c 



{bka,b\, a ) = {<WL'} = l67r 3 k+S 3 (k-k')S a 



(2.41) 



with all other commutators or anti-commutators equal zero as dictated by the spin and statistics assignments 
of Yukawa theory QED or QCD. 

The initial mass parameters mi and [J,^ do not include effects of any interactions and are independent 
of the interaction strength. We may have to consider limits where the mass parameters are close to zero, in 
comparison to all other quantities of relevance to physics. For example, /Ze°^ may be the infrared regulator 
mass denoted by [is. Recall that the subscript e indicates that the mass parameters stand in the hamiltonian 
with A = oo. 

The initial hamiltonian contains an interaction part, Q2 = Hi ~ ■ For example, electrons may 
emit photons. One writes the corresponding interaction term in QED as 

E / N^Pi]^ 3 ^ + fc - t[)v. Mi e V^CA^X- ( 2 - 42 ) 



We use conventions to be specified shortly. The hamiltonian term (2.42) is contained in the expression 



h = J dx d 2 x ± eip m (o> (x) fi(x)ip m {a) (x) 
where the fields ip (o)(x) and A v (x) for x + = are defined by writing 

tn^^) = E / W \ U m ^ka b ^ e ' lkX +V-V>ujLe ik * 
a •* 

and 



x + =0 



m ( ">kakcr L 



(2.43) 



(2.44) 



A"(x) 



E 



£k<r a k*e 



-ikx 



~ v * /it o ikx 
'-k>7 a kcr L 



(2.45) 
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Spinors u m ka and v m fc<r are defined by boosting spinors for fermions at rest, u ma and v m(T , to the momentum 
k, as if the fermion mass were m. This is done using the light-front kinematical boost representation for 
fermions 

S(m, k) = {mk + )- 1/2 [A + k+ + A_(m + a^ 1 )]. (2.46) 

Namely, u m ka- — S{m, k)u ma and v m ka = S(m, k)v ma . Solving constraint equations for the free fermion fields 
in canonical field theory amounts to using these spinors. The same boost operation defines the polarization 
vectors for photons which are independent of the photon mass. We have e\. a = (s~^ a = 0, e^ a = 2fc~ L e^:/fc + , gjL 
The spin label a denotes the spin projection on the z-axis. We adopt a number of conventions from Ref. B2 |. 
It is useful to work with the above spinors and polarization vectors because they provide insight into the 
physical interpretation of the calculated matrix elements. For example, the spinors and polarization vectors 
help in tracing cancelations which result from the current conservation (e.g. see Eq. (3.103) etc. in the next 
Section). 

Equation (2.43) includes 5 terms in addition to (2.42). The other 5 terms lead to emission of photons 
by positrons, absorption of photons by electrons or positrons, or to transitions between electron-positron 
pairs and photons. There is no term leading to creation of an electron-positron pair and a photon, or 
to annihilation of such three particles. This is the distinguished property of the light-front hamiltonians: 
conservation of momentum k + > excludes a possibility that the three momenta sum up to zero. 

Strictly speaking, one has to limit each k + from below by a nonzero positive lower bound in order to 
make sure that the three +- momentum components cannot add up to zero. This lower bound is provided 
by the inverse of the quantization volume. Our cutoffs and scale hierarchy ensure that this largest of 
spatial scales in the theory does not need to be invoked in the description of physical phenomena. The 
regularization procedure cuts off such small momenta long before they have a chance to become relevant. 
If high-order perturbation theory subsequently leads to effective hamiltonians which describe universal low 
momentum components in all physical states the notion of a nontrivial vacuum has to be taken seriously into 
account for practical computational reasons. A priori, we cannot exclude this will happen. But we postpone 
considerations of such a situation until it becomes necessary in the future work. 

The product tp finp denotes a sum of 6 basic interactions. The products of creation and annihilation 
operators are ordered as indicated at the beginning of this Section. However, Eq. (2.43) requires additional 
steps before one can assign it a well defined meaning because operators such as (2.42) can easily produce 
states of infinite norm. One needs to define the individual terms such as (2.42) in order to provide meaning 
to the whole combination of similar terms in Eq. (2.43) 

There are inverse powers of k + in Eq. (2.42) and k + may be arbitrarily close to 0. For example, when 
k% and ki in (2.42) are similar (and they are allowed to be arbitrarily close to each other no matter what 
their own size is), the photon momentum — — is arbitrarily close to zero. The problem is that 
the photon momentum appears in the photon polarization vector in the denominator: e;7 9(T2 = 2fc^~e^~ 2 /fc^~ '. 
Unless k£~ is close to zero the resulting emission strength approaches oo for k£ — > 0. Therefore, even for a very 
small coupling constant e, the interaction can be arbitrarily strong. This divergence is canceled in special 
circumstances. For example, in the tree diagrams for the S-matrix elements in QED, the cancelation is a 
consequence of the presence of more terms in the hamiltonian and the energy and charge current conservation 
in physical processes. However, for the off-energy-shell matrix elements of the T-matrix, in loop diagrams, 
or in bound state equations, such cancelations will not be ensured automatically and could lead to ill-defined 
expressions. 

In particular, one has to keep in mind that in the perturbative calculation of the S-matrix it is possible 
to apply energy and momentum conservation laws for incoming and outgoing particles on their mass-shells. 
In contrast, in the bound state calculations, the individual particle momenta cannot simultaneously be on 
the individual mass-shells and still sum up to the bound state momentum - the bound state dynamics is 
always off-shell and the on-shell perturbative mechanisms for cancelations cease to be sufficient. 

In Eq. (2.42), the inverse powers of the longitudinal momentum also appear in the fermion spinors. 
These can be a source of divergences too. However, the examples we describe in this article do not lead to 
problems with infrared fermion divergences and we do not dwell on this subject here. 

The spinor matrix elements depend on the transverse momenta of the fermions and the boson polarization 
vector depends on the transverse momentum of the boson. The strength of the interaction grows when the 
relative transverse momenta grow and leads to divergences. The divergence problem manifests itself clearly 
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when one attempts to evaluate the ratio of norms of the states h\ka) and \ka). This ratio is certainly not 
finite and it is not well defined. 

One might ask if it is useful to consider the ill-defined hamiltonian term (2.43). The answer is unam- 
biguous yes because scattering amplitudes calculated using this term in combination with two other terms 
in second order perturbation theory agree very well with observable scattering of electrons and photons. No 
loop integration appears in these calculations to indicate the divergence problem. 

It is well known that the terms one should put into the light-front hamiltonian are provided by the formal 
lagrangian density for electrodynamics C = — jF^F^ + -^{i p — m)ip. One can rewrite the lagrangian 
density into a corresponding light-front hamiltonian density by using an expression for the energy-momentum 
tensor density . Integrating T H over the light-front hyper-plane gives the expression one starts from in 
building the light-front hamiltonian for QED. 

The initial hamiltonian for QED results from formal operations on fields V>+ and A 1 - . One uses 
the gauge A + = and solves the constraint equations, substitutes expansions of the form (2.44) and (2.45) 
into the formal expression for T H , integrates the density over the light-front hyper-plane and normal-orders 
all terms. The normal-ordering produces terms that involve numerically divergent momentum integrals. The 
classical field theory does not tell us what to do with the divergences resulting from the ordering of operators. 

To deal with the divergences one has to regularize the hamiltonian theory from the outset. The naive 
connection between the classical theory and the quantum theory as given by the quantization rules is broken 
by the regularization. The regularization turns out to force new terms in the hamiltonian. To gain control 
on the regularization effects one has to construct a renormalization theory for hamiltonians. The diverging 
terms which result from normal ordering can be safely dropped in the form they appear ill-defined in the 
canonical approach because the renormalization procedure introduces other terms of the same operator 
structure to replace them. 

The regularization for light- front hamiltonians which we apply to expressions resulting from field theoretic 
lagrangian densities is first described for the term (2.42). In that term, the parent momentum P equals 
k[. The spinors and polarization vectors conveniently group a number of terms with different momentum 
dependences into a combination which is invariant under light-front kinematical symmetry transformations. 
Among those terms there are terms containing masses, terms which in field theory result from derivatives 
id or id + or from inverting the operator id + . All those derivatives are replaced in the term (2.42) by 
momenta of particles created or destroyed by that term. 

We first introduce the daughter momentum variables for the created electron and photon. We have 
introduced daughter momenta in a similar configuration in Section 2. a while defining the similarity functions 
f\. Here, we use the daughter momentum variables for the purpose of regularization. The variables are 

xi=kt/k't = x, (2.47.a) 

x 2 = kt/k't = 1 -x, (2.47.6) 

x\ = k'i/k'f = 1, (2.47.c) 

= k^ -xiP 1 - = k ± , (2.47.d) 

K^ = k^- x 2 P^ = -k ± , (2.47.e) 

n'f = k'f - x' x P^ = 0. (2.47./) 

For each creation and annihilation operator in the interaction term (2.42) we define a daughter energy 
variable. Namely, 

K l +™e « -+- 771 e 

d = = , (2. 48. a) 

Xi X 
K ±2, (0)2 12, (0)2 

X2 1 — X 
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For each creation and annihilation operator in the interaction term (2.42) we introduce a factor which is a 
function, r(yi), of the variable y\ = ee^/A 2 , where the subscript i denotes the operator in question. In the 
no cutoff limit, e — > 0. A is an arbitrary constant with dimension of a mass (Ti = c = 1). All masses and 
momenta are measured in units of A. In this article, we choose r(y) = (1 + y)^ 1 ■ Thus, the term (2.42) is 
regulated by the factor 

(1 + eei/A 2 )- 1 ^ + ee 2 /A 2 )- 1 (l + eea/A 2 )" 1 (2.49) 

under the integral. The third factor in the above expression can be replaced by 1, since is a finite 
constant and it cannot compensate the smallness of e. We shall make such replacements wherever the parent 
momentum is carried by a single creation or annihilation operator. 

In the case of terms which contain only 1 creation and 1 annihilation operator, i.e. in , no reg- 
ularization is introduced. Restrictions on the particle momenta in these terms would violate kinematical 
symmetries of the light-front hamiltonian dynamics because momenta in these terms are equal to the parent 
momenta and limiting the parent momenta violates the light-front boost invariance. 

In the initial expressions for hamiltonian densities of Yukawa theory, QED or QCD, only terms with 
products of up to four fields appear. Therefore, we have only two more situations to consider in addition to 
the cases such as and terms of the type (2.42). In the first situation we have three creation operators 
and one annihilation operator or vice versa, and in the second situation we have two creation and two 
annihilation operators. Both cases are regularized using the same general rule. 

Independently of the number of creation and annihilation operators in a product, the regularization is 
introduced by multiplying every creation and annihilation operator in the product by a function r(y) such 
as in the factor (2.49), where y = ee^/A 2 and is the corresponding daughter energy variable. Later, after 
counterterms are calculated, the same regularization factors are introduced in the counterterms. 

An additional step is required in the case of hamiltonian terms which originate from the products of 
four fields including inverse powers of id + acting on a product of two fields. We introduce two kinds of 
a fifth daughter momentum and two corresponding daughter energy variables, esi2 and 6534. The num- 
bering originates from assigning numbers to the fields in the product according to the schematic notation 
(j>i4>2(id + )~ n (j)34 l 4- One of the fifth daughter energy variables is associated with the operators coming from 
the fields number 1 and 2, and the other one is associated with the operators coming from the fields number 
3 and 4. The regularized terms will contain an additional product of functions r(y§i2) and r(j/534) with the 
arguments y 5 i 2 = ee 5 i 2 /A 2 and y 534 = ee 534 /A 2 . 

The auxiliary daughter energy variables e§i2 and 6534 are calculated as if they represented daughter 
energy variables for an intermediate particle, a boson or a fermion, created and annihilated in the vertices 
which contained the products 4>i4>2 & n d ^3^4, respectively. Those vertices are treated as if each of them 
contained three fields instead of two but the field of the intermediate particle was contracted so that the 
corresponding creation operator and the corresponding annihilation operator are absent in the resulting 
term. This particular definition of a gedanken intermediate particle does not refer to any particular Fock 
state and remains valid when the operators are replaced with q\ by the unitary transformation U\. The 
definition was inspired by Refs. [p3| and [^4| where the correspondence between the intermediate states with 
backward moving particles with spin in the infinite momentum frame and the light-front seagull interaction 
terms is extensively described. 

Mathematically, the definition of 6512 and 6534 is introduced in the following way. Every creation and 
annihilation operator in the fields </>i , fa, <j>3 and ^4 is assigned a corresponding number Sj , i=l,2,3,4. equals 
+1 for a creation operator and Sj equals —1 for an annihilation operator. We define kt = ^fc^ + s±k± \ 
and S5 = (— s 3 k^ — Sikt)/kt. The gedanken particle is thought to be created in the product of fields 
including 4>3<j)4 when S5 = 1 and it is thought to be annihilated in that product when S5 = — 1. We define 
the momentum £5 = (k^,k^) by the relation s^ks — —53/23 — s^k^ = siki + S2fc2- We also introduce two 
auxiliary parent momenta, P34 = | (fcs + ks + k^) and P12 = \ (fes + k\ + k 2 ) . Then, we introduce the daughter 
momentum and energy variables 

Z512 = k+/P+, (2.50.O) 
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e512 = ) (250c) 

2=512 

2=534 = (2.50.d) 



(0)2 



K 534 = k 5 ~ ^534-P 3 "4' (2.50. e) 



(0)2 



e 534 = K534 + m£5 , (2.50./) 

2=534 

where rnf^ equals for regularization of the terms involving {id + )~ 1 and rnf§ equals ni ^ for regular- 
ization of the terms involving (id + )~ 2 . This step completes our definition of the ultraviolet regularization 
of initial hamiltonians. 

We proceed to the definition of the infrared regularization. Inverse powers of id + for massive particles 
are already regulated when the ultraviolet regularization is imposed. This was explained above. 

For each creation and annihilation operator of an initially massless particle we introduce a factor which 
limits the daughter momentum fraction x for that operator to be greater than 5. An example of such a 
factor is given by (1 + <5/x) -1 . Note that our definition also implies that the same regularization factor is 
inserted for the gedanken particles with 2:512 defined in Eq. (2. 50. a) and X534 defined in Eq. (2.50.d). 

Besides introducing the cutoff <5 on the momentum fractions carried by massless particles, we can also 
introduce for each initially massless particle a finite regularization mass term which is denoted by m$. In 
other words, in the case of the initially massless particles, — m$. Such finite masses in the daughter 
energies lead to additional suppression of the infrared longitudinal momentum region. The additional mass 
terms are introduced through mass counterterms which contain unknown finite parts. Since the finite parts 
are arbitrary and not known to be zero we introduce the finite mass terms and investigate their role. 

2.c Renormalization conditions 

The free finite parts of counterterms are determined by renormalization conditions which result from 
comparison of theoretical predictions with data. Calculations of observables require solutions to bound state 
or scattering problems using renormalized hamiltonians. In principle, one could work with hamiltonians of 
any width A. In practice, one is limited to consider some subspaces in the Fock space. Therefore, the issue 
of setting renormalization conditions is subtle. 

In theories with small coupling constants and without confinement one has an option of defining on- mass- 
shell renormalization conditions for single particles and scattering states in perturbation theory. It means 
that one can determine free parameters in the effective hamiltonians by demanding that single particle 
eigenstates of an effective hamiltonian and the S'-matrix calculated using this hamiltonian have the required 
properties. The key examples to be discussed in detail in the next Section are Yukawa theory (pseudoscalar 
coupling) and QED. 

In theories with confinement one has to choose mass parameters for confined particles and these are 
not directly observable. We suggest in this case to use similar renormalization conditions in perturbative 
calculations of effective hamiltonians as in QED. Details are described in Section 3.c. Besides ultraviolet 
the perturbative self-interactions of quarks and gluons diverge also in the infrared region where intermediate 
states have similar energies to the outer states and the effective dynamics is no longer perturbative. Therefore, 
the perturbatively renormalized mass terms introduce large infrared effects in the effective dynamics where 
perturbative cancelation mechanisms are no longer valid. These large infrared effects are welcome as a source 
of confining potentials suggested by Perry. Jl2"| 

The key question we have to answer in practice is how many effective particles have to be taken into 
account to solve the effective eigenvalue problem and how many can be included in a doable calculation. 
A good example of a theoretical problem one can think of is how the momentum or spin of a proton 
is shared by its constituents. The phenomenology of deep inelastic scattering of leptons and nucleons 
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suggests a considerable number of constituents even at moderate momentum transfers. If the number of 
effective constituents has to be large one may encounter ambiguities in the determination of free parts of 
counterterms because observables will be calculable only through complicated procedures. On the other 
hand, the constituent quark model suggests that the large number of constituents is not needed to explain 
main features of the spectrum of hadrons. Therefore, one can expect that many states are important in the 
large width hamiltonian dynamics but only a few effective particles appear in the small width case. The 
renormalization conditions set through the small width dynamics will use a small number of constituents 
but require non-perturbative solutions for the spectrum. 

Q\ in Eq. (2.35) contains the unknown finite parts of counterterms in Ge finitei^o)- The hamiltonian 
H\ can be used to calculate scattering amplitudes and bound state properties. The most familiar example 
of QED is largely perturbative as far as renormalization is concerned in order a. One can calculate the 
physical electron energy defined as the lowest eigenvalue of the effective hamiltonian for the eigenstates with 
electron quantum numbers. Thanks to the symmetries of the light-front dynamics the eigenvalue has the 
form {p 1 - 2 + m 2 )/p + and m e has to be equal to the physical electron mass. Note also that the effective 
mass term for the interacting photons must be different from zero (and growing with A) in order to obtain 
massless photon eigenstates. Examples of the renormalization conditions for QED are presented in the next 
Section. The same procedure in QCD is expected to lead to strong infrared effects because the non-abelian 
QCD interactions prevent the same cancelation of infrared divergences as in QED. 

Hamiltonian belongs to the algebra of Poincare generators. The Poincare algebra commutation relations 
can be studied order by order in perturbation theory to find out constraints the algebra imposes on the 
counterterms. The general structure of the similarity transformation for creation and annihilation operators 
allows extension of the hamiltonian renormalization procedure to the whole algebra. The renormalization 
group evolution is given by the same Eq. (2.5) for all generators. Renormalization of the Poincare algebra 
is not further analysed in this article. |25| 



3. EXAMPLES OF APPLICATION 

This Section describes a set of examples of lowest order calculations of renormalized effective hamiltonians 
using the scheme from Section 2. We begin by the description of generic rules for calculating the right-hand 
side of the renormalization group equation (2.29). The rules follow from the commutator structure. Then, 
we discuss examples from Yukawa theory, QED and QCD. 



3. a Evaluation of commutators 

The right-hand sides of Eqs. (2.29) and (2.38) are commutators. This implies that the interactions which 
appear in the hamiltonians H\(q\) and in the counterterms in H e are connected. This Section explains how 
this result comes about. 

The commutators can be evaluated in a number of equivalent ways but some of the ways are more 
convenient than others. Suppose we are to evaluate 

H=[A,{B} d ]. (3.1) 

A = A(X, Y) 4 4 a Vj ,B = B(V, W) U[ B =1 < IL=i a™, and C = £ 2 E{z)a\a z . The right-hand 
side of Eq. (3.1) equals 

H = A(X, Y) ft al f] a V] f[ < f] 

2—1 j—1 k—1 1=1 

B(V W) Ib Jb Ia Ja 

n < n a ^(x, y) n 4, n > ( 3 - 2 ) 

fc=l 1=1 1=1 J=l 

where E w = Yli=i E( w i) an d &v = J2k=i E( v k)- By commuting n/=i a yj m tne nrs * term through JlfeLi a t k 
one generates the contracted terms with a number of contractions ranging from 1 to the smaller of the num- 
bers Ja and Ib , and a term with rX/=i a Vj standing to the right of Ofc=i a t k ■ Then, by commuting Yl.f=i a wi 
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in the latter term through Iliii a £ j one obtains new contracted terms with the number of contractions 
ranging from 1 to the smaller of the numbers I a and Jb, and a term equal to the second term in Eq. (3.2) 
with an opposite sign which thereby is canceled out leaving only connected terms in the result for H . This 
result holds despite anti-commutation relations for fermions because interactions contain even numbers of 
fermion operators. 

After the second term in Eq. (3.2) is canceled one is left with a number of partially contracted terms 
in which annihilation operators may still stand to the left of creation operators. A number of ordering 
transpositions need to be done before a generic ordering of operators adopted in the previous Section is 
achieved. In fact, the process of commuting factors in A through factors in {B}^ in the first term on the 
right-hand side of Eq. (3.2) produced above a number of terms with creation operators moved to the right of 
annihilation operators unnecessarily. These transpositions have to be undone to recover final answers with 
the adopted ordering. Nevertheless, it is visible that disconnected terms cannot appear and the following 
rule simplifies the calculations. 

The right-hand side of Eq. (3.2) equals the sum of the contracted terms which result from A{B}^ 
by moving II /=i a vj through a v k an d, the contracted terms which result from — {B}gA by moving 

fl/lfi a wi through nf=i a xi- All other terms cancel out. 



3.b Yukawa theory 



The standard procedure from Ref. p3] leads from the lagrangian density Cy — ip(i j9 — m — g<fi)ip 



i(<9 M (/><9 p — /i 2 2 ) to the light-front hamiltonian expression of the form 



Hy 



We replace fields ip m (x) and 



dx d x 



2„-L 



-d 



_L2 



2id+ 



-ipn 



+gip m ipm4' + g i> 



2id+ 



x + =0 



(3.3) 



for 



r+ - 



by the Fourier superpositions of creation and annihilation 



operators, order the operators in all terms and drop the terms containing divergent integrals which result 
from the contractions. Then, we introduce the regularization factors. 

In the course of calculating effective hamiltonians we will also add new terms to Hy due to the presence 
of the regularization, in accord with the renormalization theory from the previous Section. For example, we 



will add a small term Sm 2 



m 2 to mr in the first term and 5fij 



\x to fi in the second term. 



We will calculate these terms below using the renormalization theory to order g 2 . 

In order to consider particles with quantum numbers of nucleons and pions one needs to include the 
isospin and replace the scalar coupling by 175. jl7| However, for the purpose of the illustration of the 
renormalization procedure to second order in the coupling <?, we do not have to introduce these explicitly. 
The additional factors merely lead to somewhat different algebra which can be traced throughout the whole 
calculation and final results including isospin and 175 can be read from the results in the Yukawa theory. In 
this Section we assume m > /i > 0. 



Meson mass squared 

The simplest example of a second order expression for a term in an effective hamiltonian in the Yukawa 
theory is provided by the meson mass squared. We first describe steps which produce this expression. The 
number of distinct steps in the procedure is 10: defining the regularized initial hamiltonian, calculation of 
the effective hamiltonian, analysis of the cutoff dependence of finite matrix elements of the calculated terms 
and extraction of the structure of the divergence, evaluation of the counterterm, isolation of the finite part, 
calculation of the effective hamiltonian knowing the structure of the counterterm, solving a physical problem 
such as an eigenvalue problem or a scattering problem using the effective hamiltonian, adjusting the finite 
part of the counterterm to match data (including adjustments for the observed symmetries), and computing 
the final expression for the effective hamiltonian with the counterterm finite part determined from the fit to 
data. 
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The simplest example is described in full detail of the 10 steps. Such extensive presentation is not 
provided in later examples where more complicated expressions would require too much space. The first 
example is discussed in such detail despite the fact that in this case it is easy to predict the answer. 

For example, one might propose the structure of the counterterm using, as is usually done, some scat- 
tering amplitude instead of the matrix elements of an effective hamiltonian. Note that one can also impose 
renormalization condition using a scattering amplitude which results from a calculation performed without 
use of the effective hamiltonian. 

However, the systematic approach from Section 2 is the only tool we have for dealing with more compli- 
cated cases of light-front hamiltonians and their eigenvalue equations. In other words, the simplest available 
case is used to present all the steps in detail because it illustrates the procedure in a familiar setting. When 
we proceed to more complicated interactions details of the calculation are discussed only where a new feature 
appears. 

Equation (2.36) implies to second order in Q^x that 



lx Glx = 



Gl2X 



Gix - Eix 



J21X 



11 



£ 

p=3 



Gipx 



Gix - Eix 



Gpix 



(3.4) 



li 



where the double-digit subscripts refer to the number of creation and annihilation operators (in that order) 
and the bracket subscript denotes the part which contributes to the rate of change of Gix with A. Eix is the 
eigenvalue of Gix which corresponds to the creation and annihilation operators indicated by the subscript 
11. The reason for that only one free energy eigenvalue appears in the denominators is that Qix of Eq. (3.4) 
is a one-body operator and quantum numbers which label creation and annihilation operators in Qix arc 
the same, including momentum. Therefore, the free energy eigenvalues are also the same: both are equal to 
Eix- Consequently, all commutators are written on the right-hand side of Eq. (3.4) in the simplified form. 
The numerator similarity factors reduce to the derivative of fx (we have chosen n = 1 in Eq. (2.19)). Terms 
with more than two intermediate particles (p > 3) are of order g 4 or higher. 

Assuming that g in Eq. (3.3) is extremely small, writing Qix as a series in powers of g and keeping only 
terms order g 2 , we obtain the following result from Eq. (3.4) for the meson free energy term. 



11 meson X 



ril A;- L2 +/i? t 



(3.5) 



A remarkable feature in this result is that no correction arises to the term k ±2 /k + which is protected by the 
kinematical symmetries; the total transverse momentum docs not appear in a boost invariant expression. 
The width dependence of fix is determined by the equation 



J\XK 



ri/ 2 (z 2 )8(x-l)W 
d\ M 2 - fi 2 



I { X 1 K )? 



(3.6) 



where M 2 — (k 2 +m 2 )/a;(l — x). m 2 and fi 2 are the original bare mass squared parameters from Eq. (3.3). 
They do not include terms order g 2 and higher because such terms would lead to higher order corrections than 
g 2 for the whole expression. The terms order g 2 and higher are treated as interactions in the perturbative 
calculation. 

In terms of graphs for the effective hamiltonian calculus, Eq. (3.6) represents the contribution of a 
fermion loop on a meson line. However, the graphs are not provided in order to avoid confusion with other 
diagrammatic techniques. 



Using Eq. (2.20) in the limit of a ^-function, f(u) = 8(u — u), one obtains 



f 2 (4) 



A 2 + 



u 



U 2 



1 - \/uo 



(3.7) 



(3.8.a) 



For example, for w = \ one nas f 2 ( z x) = ^I-^ 2 + 3/i 2 — M 2 }. Therefore, the derivative of fx with respect 
to A forces the invariant mass of the fermion-anti-fermion pair, M 2 , to be equal A 2 + 3/i 2 . The derivative 
selects the range of energies in the integral where the similarity function changes most rapidly. The regions 
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where the function approaches a constant, i.e. 1 near the diagonal and beyond the hamiltonian width, are 
strongly suppressed. The region that contributes is the edge of the diagonal proximum. Jl[] The derivative 
of fx is large and positive in this region and it approaches a (5-function in the limit of Eq. (3.8). 

In the limit of an infinitesimally small uo, as discussed below Eq. (2.20), one would substitute A 2 = 



A . Then, 



f 2 (4) = o 



A + A t 2 - M' 



(3.8.6) 



The numerator factor in the square bracket in Eq. (3.6) originates from spinors of the intermediate 
fermions, Tr(j/> m + m)(p m — m) with p 2 n = p 2 n = m 2 . The subscript m indicates that the — component 
is calculated from the mass-shell condition knowing + and _L components. + and _L components of p and 
p are constrained by the light-front spatial momentum conservation law, p + p — k, where k is the meson 
momentum. The pseudoscalar interaction with ij 5 gives the same result with an additional term +8m 2 in 
the numerator. 

According to Eq. (2.49), 



r e (x,K) = 



M 2 ( M 2 , 
^v 2_+ r A 2 " ' x ^~ x > 



(3.9) 



No infrared regularization is required in Yukawa theory with massive particles, m > and /i > 0. 

If the regularization factors in Eq. (2.49) contain ej divided by 1 — Xi in place of one obtains here 



r e (x, k) = 



1 



M 

A 2 



2\ 1 



(3.10) 



instead of Eq. (3.9). The spinor bracket on the right-hand side of Eq. (3.6) equals 8(p 3 ) 2 when one changes 
variables from x and n 1 - to p with p 1 - = k 1 and p 3 is determined from the relation M. — 2\/ m 2 + p 2 . 
Thus, with the modified regularization factor of Eq. (3.10) the spinor bracket can be replaced by |p 2 = 
| (jM 2 ~ m 2 ) and the integrand function of x and k 2 can be reduced to a function of M 2 . Such simplifications 
are helpful in a qualitative analysis of the cutoff dependence. 
In the limit of Eq. (3. 8. a) for uo = \ one obtains 



d\ 2 



a 

llfrr 



X 2 + 2^ 



dz- 



o [(l + 2/a) 2 -z 2 



(3.11) 



where a = e(A 2 + 3^ 2 )/A 2 and z Q = y/l - 4m 2 /(A 2 + 3/i 2 ). Note that for A 2 < 4m 2 - 3^ 2 the derivative of 
the effective meson mass equals zero and the mass stays at the width independent value M4m 2 -3^ 2 - ^ one 
uses Eq. (3.10) instead, the corresponding result is 



dj4 

dX 2 



a 
48^ 



1 



A 2 + 2tf 



(*o)*o(l + o) 



(3.12) 



Eqs. (3.11) and (3.12) are the same for a«l which is the limit of removing the regularization cutoff, e — > 0, 
for a fixed hamiltonian width A. In this limit one has 



dX 2 



48tt 



A 2 + 2/i 2 



4m 2 



A 2 + 3/i 2 



.3/2 



6{X 2 + 3/i 2 - 4m 2 ) 



(3.13) 



If one assumes that the meson mass squared parameter in the effective hamiltonian has some finite value, 



Mo = MAn a * some such that A 2 , > 4m — 3/r then, the integration of Eq. (3.13) demonstrates that 



,^M^^(A 2 -A 2 ) + ^-(/i 2 -6m 2 )log^ 



2 

f^conv 



(A, Ac) 



(3.14) 



; (A, Ao) denotes the result of integrating the convergent part of the integrand, 



Mcon«(A, A ) 



48tt 



ds 



4m 2 



3/i 2 



3/2 



/i 2 — 6m 2 



(3.14.o) 



24 



/j, 2 onv (\, A o) has a finite (i.e. A-indcpcndent) limit for large A. It contains the terms which vanish for large A 
as inverse powers of A 2 . The dependence of M 2 on „(A, Ao) on m and fi is not indicated explicitly because we will 
not need it in the discussion of counterterm s. However, one should keep in mind that the mass parameters 
determine the value of A = y ^ m2 ~ 3/x 2 where the effective mass stabilizes. We simplify our notation 
assuming that the effective cutoffs are always above the point of stabilization. Below the stabilization point, 
the meson mass has the constant value which is independent of A. We will show later that the constant value 
is equal to the physical meson mass. 

The mass squared term in the effective hamiltonian with a non-negligible coupling g grows linearly with 
A 2 . A logarithmic correction appears with an opposite sign due to the factor Zq in Eq. (3.13), as indicated 
in Eq. (3.14). However, one cannot make contact in Eq. (3.14) with the initial hamiltonian by letting A 
grow to infinity because one would obtain an ill-defined result. The factors depending on a in Eqs. (3.11) or 
(3.12) remove the infinite growth of fi\ when A — > oo. Eq. (3.12) is simpler than Eq. (3.11) and illustrates 
the same idea so we start with Eq. (3.12). 

Equation (3.12) can be integrated over A from infinity to any finite value because the factor (1 + a)~ 4 
provides convergence for A 2 > A 2 /e. Nevertheless, the integral diverges as a function of e when e — ► 0. The 
divergence appears as a single number. Therefore, the counterterm is also a number. We add /i 2 — fx 2 to 
ll 2 in the initial hamiltonian. We also write ll 2 as a series in powers of g, li 2 = li 2 + 5fi 2 + o(g 4 ) so that 
Sfi 2 ~ g 2 . Thus, 



a 
48^ 



/i- 



s + 2/i 2 



1 - 



4m 2 



s + 3/i 2 



3/2 



A 2 



(3.15) 



This is an example of Eq. (2.33) in a perturbative application to second order in powers of g in Yukawa 
theory. The counterterm 8li 2 will be calculated following the steps described below Eq. (2.33). 

The diverging part of the integrand equals a/487r [1 + {p 2 — 6m 2 )/ s] and the remaining part is convergent. 
The convergent part of the integrand has the same structure as in Li 2 onv (X, Ao) but the integral now extends 
from A 2 to infinity instead of from Aq to A 2 . In the convergent part, one can replace the regulating factor 
by 1. Simplifications occur in the limit e — > and the result of integration in Eq. (3.15) is 



2 2 i c 2 , " 



A 2 1 



-y- + (A 2 + 3 M 2 ) + (m 2 - 6m 2 ) [loge^ + f 
-/4>m,(°° ; A ) +°(.9 4 )> 



A 2 



11 



(3.16) 



where the square bracket originates from the diverging part. 

Following the procedure described below Eq. (2.33), we define the counterterm 5li 2 as the negative of 
the integral of the diverging integrand for some arbitrarily chosen A = Ao plus an unknown finite piece 
corresponding to A and denoted by 5/j 2 j inite (X ). Namely, 



5£ = 



a 
48^ 



Y l e +(P 2 - 6 ™ 2 ) l0 § \ - ( A o + 3^ 2 ) - (M 2 - 6m 2 ) (log ^| + H 

+$V 2 efinite( X 0) + o(g 4 ), 



(3.17) 



where 



x 2 a 



/,™te( A o) - 4g7r 



A 2 

-r-ci + (m 2 



6m 2 ) C2 + (A 2 +3 Ai 2 ) + ( M 2 -6m 2 )(log^| + ^) + o(g 4 ) (3.17.o) 



with the unknown finite numbers C\ and c 2 . So, in fact, 



X 2 a 



A 2 (i + ci) + (m 2 - 6m 2 )(log i + c 2 ) 



(3.17.6) 

Since the whole expression on the right-hand side of Eq. (3.17) is merely a number, it is not necessary 
to find c\ and c 2 or any other part of it separately. One can easily find the number S/i 2 j inite (Xo) from the 
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knowledge of /i 2 at some value of A. The resulting counterterm will render well denned finite boson mass 
squared parameter in the effective hamiltonians in the limit e — > 0. Using Eqs. (3.16) and (3.17) one obtains 



A 2 

^ = p 2 + 5£ finite (X ) + ^ A 2 - A§ + (/i 2 - 6m 2 ) log ^ 



-/^(oo.AJ + ofo 4 ). (3.18) 



Equation (3.18) is an example of Eq. (2.35). 

The unknown finite term 5(i 2 f inite (\o) has to be found by comparison with data, which might include 
symmetry requirements. We shall discuss an example of a renormalization condition later in this Section. 

Let us assume now that at some arbitrarily chosen value of A = Ai the effective meson mass squared 
required in Eq. (3.5) by a fit to data equals /i 2 , i.e. = p\. We can calculate Sfi 2 j inite (\o) using Eq. 
(3.18) with A = Ai and p 2 Xi on the left-hand side replaced by the number \j\ inferred from the experimental 
data. The result is 



c- 2 l \ \ 2 2 ^ 



A 2 -A 2 + ( M 2 -6m 2 )log^ 



+ £ onv (^,\ 1 )+o(g i ). (3.19) 



Note that one has to include the contribution of the convergent terms in the determination of the arbitrary 
constants when using the value of /i 2 . Knowing 5/j, 2 f inite {M) one can calculate /i 2 . Namely 



,2 _ ,2 



a 



= Ml + 48. 



A 2 - A 2 + (^ -6m 2 ) log ^ 



» 2 conv (\,\i) + o(g 4 ). (3.20) 



As expected, Eq. (3.20) is the same as Eq. (3.14) when \x = A and = fi . One can also trace the origin 
of all the terms from Eq. (3.14); the diverging and converging terms and the counterterm in Eq. (3.15). 

We can now analyze Eq. (3.11) analogously to Eq. (3.12) without calculating all integrals explicitly. 
Integrating both sides of Eq. (3.11), we have 

poo i 2 

A=fi- dsfL. (3.21) 
J\2 as 

By demanding that /i 2 removes the diverging (i.e. e-dependent) part of the integral, and leaving the finite 
part of it 2 free so that at some A = A the effective boson mass squared parameter has a desired value /Xq , 
we obtain 



L 



,4 - //,* + I ds^. ,3.22) 



The integrand is given by the right-hand side of Eq. (3.11) with A 2 = s. Since s ranges only from A 2 , to A 2 
and both are finite we can take the limit e — > under the integral sign and the integrand becomes equal 
to the right-hand side of Eq. (3.13) with A 2 replaced with s. Integration over s produces exactly the same 
answer as given by Eq. (3.14). Thus, we see that the result of Eq. (3.14) is independent of the regularization 
scheme. It is determined by the initial hamiltonian of Yukawa theory as given by Eq. (3.3). The only 
unknown in Eq. (3.14) is the value of /i 2 ,. More precisely, we know that /i 2 , = ji 2 + 5/1q + o(g 4 ) and the 
unknown term is (5/i 2 , <~ g 2 . 

Note that the above calculations can be carried out in a different way using the following observation. 
Equation (3.6) in the lowest order of perturbation theory has a particularly simple structure. Namely, the 
only dependence on A comes from the similarity function f{z 2 ) and both sides of the equation are equal to 
the derivatives with respect to A. Therefore, one can directly integrate both sides and obtain a compact 
integral expression for /i 2 for arbitrary functions /. 

One should remember that such simplifications do not occur in higher order calculations or beyond simple 
pcrturbative expansion when the coupling constant and mass parameters depend on A themselves. Therefore, 
we stress that the dominant contribution to the rate of change of /i 2 with A comes from the edge of the 
diagonal proximum. This fact remains generally valid and the procedure applied above represents a generic 
situation despite simplicity of the example. This example does not involve a distinction between the bare 
coupling g and a renormalized coupling because to order g 2 there is none. 

In order to determine (5/Zq we need to specify a renormalization condition. A natural requirement for 
/Xg is that the effective hamiltonian has one boson eigenstates parametrized by momenta p + and p 1 - with 
eigenvalues of the form p~ = (p ±2 + jl 2 )/p + where ft is equal to the physical boson mass. Our approach 
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preserves kinematical symmetries of the light-front frame explicitly and the eigenvalue is guaranteed to 
appear in that form. Therefore, one can calculate a whole spectrum of eigenvalues for eigenstates of different 
momenta by calculating the single mass parameter jl. 

In order to write the effective eigenvalue equation and hnd out (Mq which leads to the desired value of jl 
(if it is possible), the following steps need to be taken. 

One inserts Eq. (3.14) into Eq. (3.5). Then, one replaces the bare operators at and in the whole 
Q\ by the effective ones, a\ k and a\k, obtaining G\. Next, by applying the operation F\, one obtains the 
effective hamiltonian with the form factors in the vertices, H\ = F\\G\\. 

The effective eigenvalue equation for bosons is an example of Eq. (1.2). Here, it is sufficient to consider 
the eigenvalue equation for H\ in the expansion to second order in g since we are interested in <5/Iq which 
is proportional to g 2 . The only terms that contribute are the free energy term including the effective mass 
squared and the interaction terms which change the particle number by one. The latter equal the canonical 
interactions with the similarity form factors. 

To zeroth order in g, a physical meson state equals a single effective meson state, and jl 2 = fi 2 . 

No terms order g arise in jl 2 and the next correction is order g 2 . This correction can be calculated using 
the operation R and the model hamiltonian defined in a perturbative expansion from Eq. (1.3). Using 
expansion into a series of powers of g to second order, one can restrict the model space to the single effective 
boson sector. The effect of coupling to the fermion-anti-fermion pair states is included in perturbation theory. 

Another method is to arbitrarily limit the number of effective Fock sectors and diagonalize the effective 
hamiltonian in that limited space. Such procedure could be called the effective Tamm-Dancoff approach 
(ETD), cf. Rcfs. H and fl0|| . The term of the second order in g in the eigenvalue will determine 8[1q. 
One can limit the space of states to one effective boson and effective fermion-anti-fermion pairs. Note that 
interactions in the ETD are tempered by the similarity factors of width A which is on the order of particle 
masses and no ultraviolet divergence exists. This way our ETD approach overcomes the old problem of 
ultraviolet mass renormalization in the TD approach. 

The model calculation using R and the ETD calculation have to agree with each other for small coupling 
constants in the presence of a finite energy gap between states with different numbers of particles, which is 
the case here for < /i < 2m. We discuss only second order corrections in g. Therefore, we can focus on 
a straight-forward perturbation theory anyways. Nevertheless, our simple calculations have two interesting 
aspects. 

The first one is that no coupling renormalization effects arise to order g 2 . Therefore, the expansion in 
powers of g up to g 2 is equivalent to the expansion in powers of an effective coupling, no matter how the 
latter is defined. It is important to realize that the expansion in powers of g is understood here to be the 
expansion in powers of the effective coupling which appears in the effective hamiltonian of width A; g\ = g 
to order g 2 . It is not meant to be the expansion in the initial field theory coupling constant. 

The second aspect is following. The perturbative expansions applied in the effective eigenvalue problem 
are expansions in the interaction which is suppressed in strength by the similarity vertex form factor of width 
A. Therefore, the range in momentum space of the effective interaction terms is infinitely smaller than the 
momentum range of the analogous interaction in the bare hamiltonian. In other words, the effective strength 
of the interactions is greatly reduced and much smaller than the value of g itself would imply if it stood 
in the initial bare hamiltonian. Our initial expansion in powers of g can now be understood as a shortcut 
to the expansion in powers of the effective coupling. The latter expansion may have a considerable range 
of rapid convergence because the form factors reduce the size of coefficients in the expansion. The effective 
coefficients are much smaller than in the case of the initial hamiltonian without form factors. 

Thus, the operation R on H\ expanded in powers of the effective interaction (the coupling constant itself 
can be sizable), opens new options for studying the effective eigenvalue problem in the whole Fock space 
using the basis built with the effective creation and annihilation operators corresponding to the width A. 
One can estimate contributions of various components by making different choices of the model spaces and 
solving model dynamical problems numerically. Wave functions are expected to fall off sharply for large 
momenta and particle numbers if g is not too large. The effective coupling constant needs to be set equal to 
the right value at A. 

The second order expression in perturbation theory implies 

P±2 p t M V b) = ^t^V lP) - {p'\Fx[G 12X } — - * + F x [G 2 ix]\p) ■ (3.23) 
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\p) denotes a single effective meson state with momentum p + and p^, (p'\p) = 16-K 3 p + S 3 (p' — p). The term 
order g 2 produces 



~f=A-\ l**]9f&) %p^fi 9f(4) + o(9% (3-24) 
where the notation is the same as in Eq. (3.6). Using Eq. (3.8) with u = \ at A = A one obtains 

For A < \J Am 2 — ip? the effective meson mass parameter equals the physical meson mass, as promised. 
For fermions with masses order 0.9 GeV, this implies no corrections to a light meson mass such as fj,„ 
for cutoffs smaller than 1.8 GeV. But one has to remember that the correction for the pseudoscalar ttN 
interaction is different from Eq. (3.25), i.e. the spin factor has to be enlarged by 8m 2 (see comments above 
Eq. (3.9)). 

However, the actual measure of the off-shell effects is not given directly by $ but by the sum of [Iq and 
the self-energy resulting from the effective interactions. According to Eqs. (3.22) and (3.24), the sum of 
both contributions in the physical pion mass itself is actually equal zero to order a. 

Using Eqs. (3.21), (3.22) and (3.25), one can express the meson mass squared term in the initial renor- 
malized hamiltonian in terms of the physical meson mass ft and the initial mass parameter p. Namely, 
/j, 2 = /i 2 + o(g 2 ) and 

A = + 9 2 J [x K ] 8{X M2 ^ 2 r e (x, k) + o(g 4 ) . (3.26) 



Fermion mass squared 

In complete analogy to Eqs. (3.4) to (3.6) one obtains the fermion energy operator, 



hfermionX = E/ W ~ feT^ L^a + d kM- (3.27) 



Results for fermions and anti-fcrmions are identical. We have 

(p'm + m)u 

■ma- k 



J [x /W^+^ r t (x, «). (3.28) 



dX 

M 2 = {m 2 + n 2 )/x + (p 2 + k 2 )/(1 - x). The regularization factor of Eq. (2.49) implies 



r e (x,K) = 



i + ^m 2 + P 2 k2 + w2k2 + ^ 2 



,A 2 J 

The spin factor in Eq. (3.28) can be rewritten as 



(3.29) 



Umakip'm + m)u mak = u mak [x $ m + to + ]p + [p m - xk m )]u mak . (3.30) 

Jt m between spinors is equivalent to to. The term with j + is typical in light-front calculations. Its part 
proportional to k ±2 /k + cancels out. The term linear in k 1 - docs not contribute because it is odd in k 1 - and 
all other factors including the regularization factor depend only on the modulus of k^. The spin factor is 
thus reduced to 



(x + 1)to + 



7+ k 2 + (1 - x 2 )m 2 
2k+ x 



K 2 + (l + x) 2 m 2 

U m ak = • (3.31) 



Result for a pseudoscalar interaction with ij 5 is the same except for the opposite sign in front of x in the 
numerator. Inclusion of the isospin introduces the number of bosons in the theory in front of the integral in 
Eq. (3.28). 
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We observe a similar structure in Eq. (3.28) as in the meson mass dependence on A 2 in Eq. (3.6). 
Namely, there are terms diverging linearly and logarithmically and there is a series of convergent terms. We 
observe that the divergences amount to a number which grows when e — > and integrate both sides of Eq. 
(3.28) to obtain 

where, according to Eq. (2.18), z\ = {M 2 — m 2 )/(A4 2 + m 2 + A 2 ). The e-dependent terms originate from 
1 in the bracket which is independent of A. The counterterm Sm 2 <~ g 2 in m 2 = m 2 + Sm 2 + o(g 4 ) removes 
the divergence. The finite part of the counterterm is left to be determined by data. 

ml = m 2 + Sm 2 tfmlte+ g 2 J [ XK ]f(zD ^f^ff ' ^ 

The value of Sm 2 j inite is determined from the value of m 2 required in the effective hamiltonian H\ by some 
physical condition. If we had defined the divergent part by an integral from Ao we would have to take into 
account that m 2 j inite depends on Ao to compensate for the Ao dependence of the diverging integral. When 
we define the counterterm to be given by the whole A-independent part of the integral in Eq. (3.32), plus 
a finite constant to be determined by data, then 5m 2 ^ inite does not depend on Ao. Nevertheless, it can be 
expressed in terms of m\ a . For example, if the effective fermion mass squared at A = Ao should be m\ a = m 2 , 
then 

ml = m 2 + . 9 2 /m [f 2 (z 2 ) - f 2 (z 2 o ] K [ + { ^ 2 + Sf - ( 3 - 34 ) 

to 2 , = m 2 + 5ml + °(ff 4 )- ^ m o ~ S 2 can be found from a renormalization condition for the physical fermion 
mass. 

A natural condition for fitting TOq is that the effective hamiltonian at the scale Ao has fermionic eigenstates 
with eigenvalues of the form p~ = (p ±2 + m 2 )/p + , where m denotes the physical fermion mass. In analogy 
to Eqs. (3.23) and (3.24) one obtains 

~2_ 2 ft i tr 2^ 2 + i 1 + x) 2 m 2 2 



rh z =<-J [™\gf(4) x{M 2_ m 2 ) am) + o(g 4 )- (3.35) 

So, 

ml = - 2 + g 2 J M/ 2 (^J ^ 2 +y + o(g% (3.36) 

The initial m 2 can be calculated in terms of m 2 , m 2 , g and e from Eq. (3.32). The effective fermion mass 
parameter in the interacting hamiltonian of width A is 

ml ^ m 2 + g 2 J [ X n]f\zD + o(g% (3.37) 

Analogous equation in the case of nuclcons coupled to pions is the same as Eq. (3.37) except for (1 — x) 2 
instead of (1 + x) 2 in the numerator and the isospin factor 3 in front of the integral. In the limit of a 
^-function for the similarity factor / one obtains the result m 2 = m 2 + o(g 2 ), rh — mjv, and 

m 2 = m 2 N + 3g 2 J [xn]0(X 2 + 3m 2 - M 2 )^^f + o(g% (3.38) 
For A 2 = (m + n w /i^) 2 — 3m 2 , where is a small integer one obtains (a = .g 2 /47r) 

m l = m N + m %T~ c ■ (3.39) 

47T 

c = 4/3(n7r /iTr/mjv) 3 + o(fi^). The expansion formula for c shows the correction is small for small meson 
masses. Note that A 2 must be negative for small n v , as explained below Eq. (2.20). The exact result for 
n v = 3 gives c ~ 0.03 and — 4 gives c <~ 0.12. Even for quite large couplings the effective mass parameter 
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in the hamiltonian deviates only a little from the physical nucleon mass if the hamiltonian width allows 
momentum changes of the order of a few meson masses only. Moreover, the physically relevant off-shell 
effects are not given directly by the above numbers but by the difference between these and the effects of the 
interactions present in the effective hamiltonian. The combined effect is zero for the nucleon mass itself to 
order g 2 . Eq. (3.39) suggests that the self-interaction effects can be calculable in perturbation theory. This 
is encouraging for the program outlined in Ref . ||l7f . 

If we used Eq. (2.20) in the ^-function limit with an infinitesimal uq and A 2 = X 2 / ^/uq for n = 1 the 
^-function under the integral in Eq. (3.38) would be replaced by 8(X 2 + m 2 — M 2 ). 



Fermion-fermion interaction 

Our next example is the second order calculation of the effective hamiltonian term which contains prod- 
ucts of two creation and two annihilation operators for fermions. The differential equation we need to 
consider results from Eq. (2.29) for the two-fermion terms; 



d G 
a\ 



(3.40) 

J 22 

The subscript 22 denotes a term with two creation and two annihilation operators for fermions. 21 denotes 
a term with one annihilation operator and one creation operator for fermions and one creation operator for 
bosons. 12 denotes a term which annihilates a fermion and a boson and creates a fermion. For a hcrmitean 
hamiltonian, we have Q12 = Q\\- 

The right-hand side of Eq. (3.40) does not contain disconnected interactions (it never does, cf. Section 3. a) 
and one can isolate the terms with two creation and two annihilation operators for fermions by contracting 
one creation operator and one annihilation operator for bosons. The only term which contributes is the 
ordered and regularized third term on the right-hand side of Eq. (3.3). Thus, in Eq. (3.40), we have 

Q2W = [fcifcfc 2 ]167r 3 (5(fci + k - k 2 ) gu m k 1 a s u m k 2 ar{eef / K 2 )r{eeb/ K 2 )b\ iai a\bk 2 a 

= J [ p \j^ J [ XK \9u m xP+ K a f u m Pa r(ee f /A 2 ) r(ee b /A 2 ) fei P+Ka/ a[ 1 _ x)p _ K 6 Pcr , (3.41) 



no- 



where e/ — (k 2 +m 2 )/x and e b = (k 2 + /i 2 )/(l — x). This representation illustrates appearance of the parent 
and daughter momentum variables in the interaction term. The factor g u m k 1 a f u m k 2 cr r ( ee / f^ 2 ) r ( ee b/-^ 2 )bk iaf 
will be denoted by 921- The analogous factor in Q12X will be denoted by g\2- Similarly, in the case of the 
four-fermion interaction term, we have 

&2A = /t P ]p+ H /M[w]522A6ip +KCTl &J 1 _ x) p_ KCT3 &s,P+pa 2 &(l-j / )P-pa 4! (3-42) 

CT1(T2(T3(T4 

where 322 a is a function of the daughter momentum variables x, k^, y, and the fermion spin projections 
on the z-axis: o\, 02, 03 and 174. Details of the notation will become clear shortly. To order g 2 , only f\ 
depends on A on the right-hand side of Eq. (3.40) and it can be written in terms of the coefficient functions 
as 

^322 = [/{-/'} - {-/'}/] [.9l2.92l] 2 2 ■ (3-43) 

The subscript A and arguments of the functions are not indicated, to simplify the formula. Expression in the 
first bracket is called the inner similarity factor for Q22 ■ The word "inner" is used to distinguish this factor 
from the form factor introduced by the operation F\ when this operation is applied to G22 ■ The latter form 
factor can be called the outer similarity factor because it depends on the incoming and outgoing invariant 
masses only, independently of the internal structure of the interaction. 

Momentum variables in Eq. (3.43) can be expressed by the daughter momentum variables from Eq. 
(3.42). One needs to express the parent and daughter variables of G12 and Q21, and the energy denominators, 
in terms of x, K , y and from Eq. (3.42). This is done as follows. 
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The momentum labels of the fermion annihilation operators in Eq. (3.42) arc denoted by k 2 and k\ 
and the momentum labels of the fermion creation operators in Eq. (3.42) are denoted by fci and kz . The 
numbers assigned to the fermion momenta correspond to the numbers labeling their spin projections on 
the z-axis in Eq. (3.42). In terms of the origin of the annihilation and creation operators in Eq. (3.42), 
Qi\ provides the fermion creation operator with momentum k\ and the fermion annihilation operator with 
momentum k 2 . G12 provides the fermion creation operator with momentum k 3 and the fermion annihilation 
operator with momentum k 4 . There is a change of sign due to the reordering of the fermion operators. The 
boson operators from G12 and G21 are contracted and provide factors similar to the factors associated with 
an intermediate particle in the old-fashioned hamiltonian calculations of the S-matrix. 

It is useful to think about the effective hamiltonian in terms of a scattering amplitude with two vertices 
but the reader should remember that the formula we are describing is not for an S'-matrix matrix element. 
Therefore, the "scattering" language has a limited meaning. The fine point is that, after evaluation of Gx, one 
has to replace the bare creation and annihilation operators by the effective ones in order to obtain G\. The 
term in G\ is not directly related to any scattering process before the replacement is made. The replacement 
prevents confusion between the hamiltonian calculus which uses the bare operators, and the S'-matrix calculus 
which uses the effective hamiltonian and the corresponding incoming, outgoing and intermediate states of 
effective particles. The scattering language becomes particularly confusing in theories with gauge symmetry, 
spontaneous symmetry breaking and confinement. None of these features appear here in the calculation to 
order g 2 . Therefore, the scattering language is useful in the current example. 

The intermediate boson momentum is defined for + and _L components as k§ = k 3 — k^ and k§ = 
(k^ 2 + fi 2 )/k^ . These four components form the four-momentum of the exchanged boson. It is denoted by 
k^^ to indicate the mass which determines the minus component. The same result for k^^ is obtained by 
subtracting k\ from k 2 instead of £4 from k 3 . It is so because the translational invariance of the hamiltonian 
on the light-front implies momentum conservation for the + and _L components. 

Thus, the inner similarity factor in Eq. (3.43) is 

[/{-/'} - {-/'}/] = f(z! 2X ) l ~ f lf 2 f (3-44) 

The prime denotes differentiation with respect to A. 

In the case of Gi in Eq. (3.4), the whole inner similarity factor of the analogous structure was equal 
to the derivative of f 2 (z 2 ) divided by a single denominator. For both functions / in Eq. (3.4) had the 
same argument z 2 and the two corresponding energy denominators were the same. In Eq. (3.44) we have 
two different arguments of the similarity functions / and two different energy changes. Namely, z 42 \ and 
AE i2 for the vertex of G12 with momenta k 3m , k± m and £5^, and z 2 \\ and AE 21 for the vertex of G21 with 
momenta k\ mi fem and &5 M . 

The parent momentum for the vertex of Q\2 is P\ 2 = (k^^ + k 3m + k 4m )/2 so that for the + and _L 
components we have P\ 2 = k 3 . Similarly, the parent momentum for the vertex of G21 is P21 — (^5^ + k\ m + 
^2m)/2 so that for the + and _L components we have P21 — k 2 . 

Now, the rules provided by Eqs. (2.12) to (2.20) imply the following formulae for the arguments of the 
similarity functions /. 

AM 2 12 = (k bll + k 4m ) 2 - k\ m = 2(fc 5(U + k 4m - k 3m )P 12 . (3.45) 



EAi? 2 = M\ 2 + 2m 2 . 



z\ 2 \ 



AM ? 9 



SM? 2 + A 2 ' 



AM 2 21 = k 2 m - (fc 5M + k lm ) 2 



-2(k. 



5ju ~r Mm 



^2m)^21- 



(3.46) 

(3.47) 
(3.48) 



EM , = -M01 + 2m 



Z21X ~ EM 2 



AM 2 ! 



21 



A 2 ' 



(3.49) 
(3.50) 
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Equations (2.12) to (2.16) imply 



AE 



12 



and 



AE 2 



AAi 2 2 

P+ 

^12 



P+ ' 

^21 



In terms of the familiar parameters x, k^, y and p 1 - from Eq. (3.42), i.e. 

P = ki + k 3 = k 2 + k 4 , 

x = k+/P+, 



and 



kf - xP A 



y = k+/P + , 



p 1 - = k^ - yP 1 ^ , 

the mass differences which determine the arguments of the similarity functions read as follows. 



AMj 2 = {l-x) 



(k 1 - - p^) 2 + p 2 p ±2 + m 2 K ±A +m 



,2 1 



+ 



AM 2 21 = y 



p ±2 + in 2 

y 



y — x 1 — y 1 — x 

K ±2 +m 2 ( K ± _ p ±^2 + ^2 



y-x 



The corresponding energy denominators are 



and 



AE 12 = 



AE 21 = 



(k 1 - — p^Y + M P +m k + m 



y-x 



+ 



l-x 



p ±2 + m 2 K JL2 + m 2 ( K ± _ + ^ 



!J 



y-x 



/P+. 



(3.51) 
(3.52) 

(3.53) 
(3.54) 
(3.55) 

(3.56) 
(3.57) 

(3.58) 
(3.59) 

(3.60) 

(3.61) 



In Eqs. (3.40) to (3.61) always y > x. Evaluating matrix elements of the effective interaction between 
states of indistinguishable fermions one obtains results in which the momentum and spin variables are 
properly symmetrized (antisymmetrized) as dictated by the statistics. 

Evaluation of the second bracket in Eq. (3.43) gives 



[5i252i]22 = - gu m (i-x)p-Ka 3 Um(i-y)P-pa 4 r(ee 4 /A 2 )r(ee 12 /A 2 ) 
i 

■ gu mxP+Kai u my p +pa2 r(eei/A 2 )r(ee 2 i/A 2 ) . 



{y - x)p+ 

The arguments of the regularization factors appear in the mass differences. Namely, 



AMj 2 =e 4 + e 12 - m 2 = K ™ + — + + ^ 



where 



Xl2 



Xl2 



i-y 

l-x' 



1 - £12 



m 



(3.62) 



(3.63.a) 



(3.63.6) 
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Ki 2 = -p 1 - + X i2 K A 



and 



where 



AM 



to 



ei + e2i 



m 



,±2 
l 21 



X21 



k ±2 4- II 2 
1 - X 2 1 



«21 = -, 

2/ 



v 21 



a^2iP 



_l 



(3.63. c) 



(3.64.a) 



(3.64.6) 



(3.64.c) 



Eqs. (3.62) to (3.64) contain daughter energies in the notation introduced already in Eqs. (2.48) and (2.49). 
Similar subscript notation will be used to label the regularization factors. Equation (3.43) reads, in the 
abbreviated notation, as follows. 



dX 



922 



V /12/21 _ 1 - x /{2/21 
- x AMn 



y — x AM\ 2 



9{y - x) gu 3 u 4 r 4 r 12 gu x u 2 nr 2 i 



(3.65) 



In the familiar limit of Eq. (3.8) where the similarity function / approaches the ^-function, i.e. /12 = 
8[X 2 + 2m 2 - AM 2 2 \ an d /21 = 9[X 2 + 2m 2 + AM 21 ], the derivatives of the similarity functions become 
(5-functions and one can integrate Eq. (3.65) using the relation ds9(s - b)6(s - c) = 9(c - b)6(c - a). The 
result is 



522A = 522e + 



2/#21-12#21-A (1 - X)9 12 -\0i2-21 



_(y-x)\AM 2 21 \ (y-x)\AM 2 2 
The symbols 9 with various subscripts denote the following functions: 

^12-21 = 1 - 021-12 = 0(\AM\ 2 \ - \AM 2 2l \), 



9(y - x) gu 3 u 4 r 4 r 12 gu x u 2 r x r 21 . 



(3.66.a) 



(3.66.6) 



9 12 ^=9(\AM 2 12 \-\ 2 -2m 2 ), 



and 



(3.66.c) 



721-A 



9{\AM 2 2l \-\ 2 -2m 2 ). 



(3.66.d) 



The initial value term of g 22e at A = 00 is absent in the canonical hamiltonian. It is equal zero if matrix 
elements of H 22 \ = F\ [G22A] between finite free energy states have a limit when e — > 0. If the limit does 
not exist due to a diverging e-dependence a counterterm containing nonzero g 22e is required to remove the 
divergence. 

If we used Eq. (2.20) in the ^-function limit with an infinitesimal u and A 2 = X 2 /^/uo for n = 1 then, 
A 2 + 2to 2 in the 9- function arguments above would be replaced by A 2 . This feature will be used later in the 
case of nonrelativistic bound states. 

The easiest momentum configuration to analyze g 22ti is the one where the sum of free energies for the 
momenta of creation operators equals the sum of free energies for the momenta of annihilation operators: 
(ki m + k 3m ) 2 = (k 2rn + k 4m ) 2 = M 2 . We will refer to this configuration as the energy- diagonal part of the 
interaction. In the energy-diagonal part of the interaction, we have AE\ 2 = —AE 2 \, and 



,12 



m 



,±2 



m 



,±2 



m 



n-L2 



TO 



i-y 



l-x 



Thus, 



\AM 2 „ 



\AM 2 21 \ 



V 2 + q 2 

y-x 



(3.67) 



(3.68) 
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where q 2 = (k 1 - — p^-) 2 + (y — x) 2 M 2 . These relations imply 6*12-21 = 1 — #21-12 = 9(1 — x — y), 612-x = 
9[p 2 + q 2 -(y-x)(X 2 + 2m 2 )/{l-x)} and 621-x = 9[p 2 + q 2 - (y - x){\ 2 + 2m 2 )/y}. Therefore, 6 12 -x = 1 in 
the same momentum range where 621- x = 1- Thus, the energy-diagonal part of the fermion-fermion effective 
interaction order g 2 is 



gu z u^r^rx2guxU2r x r2\ . 

322 A = ff22e H 2 , ^ 2 % - z)0 

M + 9 



2 1 -2 

m +q 



max(y, 1 — x) 



(A 2 + 2m 2 



(3.69) 



Equation (3.69) is helpful because it provides insight into the more complicated interaction from Eq. 
(3.66). When the momentum transfer is sufficiently large and (y — x)(\ 2 + 2m 2 )/max(y 1 1 — x) is negligible so 
that the f-functions and the regularization functions n, 7*12, r\ and r 2 i in Eq. (3.69) equal 1, the second term 
on the right-hand side of Eq. (3.69) is equal to Feynman's expression for the one boson-exchange scattering 
amplitude for two fcrmions. Namely, the numerator factors are standard for the Yukawa interaction and the 
denominator equals p 2 — [k^ m — fct m ) = p 2 — (&2m — ki m ) 2 (the necessary antisymmetrization for identical 
fermions requires evaluation of matrix elements of the hamiltonian term under consideration). However, 
there is a difference between the energy diagonal part of the effective hamiltonian matrix elements and the 
on-shell Feynman scattering amplitude due to the ^-functions and the regularization factors in Eq. (3.69) 
where they differ from 1. In other words, our theory is regularized ab initio and the resulting amplitudes 
contain the regularizing factors. The width dependent factor of the effective hamiltonian does not belong in 
a physical scattering amplitude and we will explain how it goes away when one evaluates S'-matrix elements 
later in this Section. 

The 9- functions force the momentum transfer carried by the intermediate boson, \q\, to be larger than 

1 /2 

[(y-x)(X 2 + 2m 2 )/y- p 2 ] 1 . The size of this number depends on how large A 2 and the ratios of the 
longitudinal momenta are. If A 2 is negative and compensates 2m 2 , the lower bound on the momentum 
transfer is absent. For larger A 2 , the ratio of y — x to the parent fermion y has to be smaller than p 2 (X 2 + 
2m 2 )~ 1 for the lower bound on the momentum transfer to be absent. Otherwise, the momentum transfer is 
limited from below. This means that the effective interaction term does not include the long distance part 
of the Yukawa potential. 

The regularization factors in the limit e — > converge pointwise to 1. No diverging cutoff dependence 
is obtained when evaluating matrix elements of G22A between states of finite invariant masses M 2 so the 
matrix elements of H22X arc free from divergences. Therefore, </22e = 0. We can replace the regularization 
factors in the limit e — > by 1 . 

We proceed to the analysis of Eq. (3.66). No divergences appear in the finite matrix elements of the 
effective hamiltonian of width A when e — > 0. One can see that this is the case using Eqs. (3. 63. a) and 
(3. 64. a). Namely, the arguments of the regularization factors are finite for finite A.M 2 2 and A.M21 and 
they approach when e — > 0. One demands that the free invariant masses of the states of fcrmions used to 
calculate the matrix elements are finite. The only possibility for A_M 2 2 or AA^i to diverge emerges when x 
approaches y, i.e. when the longitudinal momentum transfer between the fermions approaches zero. In such 
case, ei2 and e2i approach infinity even for a vanishing transverse momentum transfer because the meson 
mass squared is greater than zero. 

Now, the remaining factors of spinors and energy denominators, the latter multiplied by the boson phase- 
space factor of y — x, are all finite in the limit x — > y. The regularization factors r i2 and r 2 i deviate from 
1 only in the small region in the momentum space where \x — y\ < eh 2 / p 2 (or in a still smaller region 
for a nonzero meson transverse momentum). All other factors in the interaction are finite in this region. 
Therefore, for finite wave packets or bound state wave functions used in the evaluation of the matrix element, 
this small region produces a contribution which is proportional to e. Thus, it vanishes in the limit e — > 0. 
Consequently, the matrix elements of <?22e are equal and the regularization factors can be replaced by 1. 

The full result for the effective fermion-fermion interaction in the limit e — > is 



H22X = F\ [G 2 2a] 



[P] 



P+ 



J M lyP\ 922xf{zl 2 x)b\ 



XxP+K,a 1 X(l-x)P-Ka 3 0x V p +P a 2°X(l-y)P- P a 4 , 



(3.70.a) 



where 
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922\ 



712-2ltf21-A P12-Af21-12 
/i 2 - /J 2 - g 2 2 



6*(y — x)guzUigu\U2 , (3.70.6) 

912 = hm - hm, (3.70.C) 

921 = k 2m - hm, (3.70.d) 
^12-21 = 1 - ^21-12 = [(1 - x)(p 2 - q\ 2 ) - y(p 2 - q 2 21 )] , (3.70.e) 
012-x = 0[{l- x)(p 2 - q\ 2 ) -(y- x)(X 2 + 2m 2 )] , (3.70./) 

021-A = ^ [y(^ 2 - 9 2 2 i) - (2/ - ^)(A 2 + 2m 2 )] . (3.70. 5 ) 

The argument of the outer similarity factor /(zf 2A ) that limits the width of the effective interaction in 
momentum space, is 

The mass difference, AA4 22 = -M 24 — -M 2 3 , and the mass sum, £.M 22 = -^24 + -A4f 3 , are expressed by the 
fermion momenta through relations -M 24 = (k 2m + k im ) 2 and A4 2 3 = {k\ m + fc 3m ) 2 . 

Equations (3. 70. a) to (3.70.h) explain the structure of the fermion-fermion effective interaction order g 2 
in terms of the two four-momentum transfers, 912 and g 2 i- The transfer g 12 appears in the vertex where the 
intermediate boson is annihilated and the transfer q 2 \ appears in the vertex where the boson is created. The 
0-functions exclusively select which momentum transfer appears in the denominator. The lower bounds on 
the momentum transfers depend on the ratio of \x — y\ to y and 1 — x and on the masses and A 2 . 

We can now evaluate matrix elements of the T-matrix between effective two-fermion states using the 
hamiltonian of width A to second order in g; 

T{E)=H IX +H IX - 1 , H IX . (3.71) 

E - iioA + i£ 

We have H ox = Gi X and H IX = F x [G 22X + G\ 2X + G 2 i\]- The first term on the right-hand side of Eq. (3.71) 
contributes solely through H 22X . In the second term, only Hi 2X + H 2 \ X contributes in Hi X . 

The first term in Eq. (3.71) has its matrix element given by the antisymmetrization of the right-hand side 
of Eq. (3.70.b). The multiplication by f(z 22X ) does not matter because f(z 22X ) = 1 in the energy-diagonal 
matrix elements and only the energy-diagonal part contributes to the cross section. The energy-diagonal 
part of g 22X is given by 

gu 3 u 4 guiu 2 

ff22A = ~ - 2 — % -x)6 x , (3.72 

[i A — q z 

where 8 X = 8 \max(y, 1 — x)(p 2 — q 2 ) — (y — a;) (A 2 + 2m 2 )] and q = q\ 2 = g 2 i- The antisymmetrization of 
the right-hand side of Eq. (3.72) produces the contribution of the first term in Eq. (3.71) to the scattering 
amplitude. 

The second term in Eq. (3.71) provides the one-boson exchange amplitude with form factors in the 
fermion-boson vertices. The form factors are the similarity functions f x . The resulting amplitude is given by 
the antisymmetrization of Eq. (3.72) with 9 X replaced by the product of the form factors. In the ^-function 
limit, the vertex form factors equal 1 — #21- a and 1 — #i2-A- Their product equals 1 — 6 X . Thus, the second 
term provides the same contribution as the first term but the factor 9 X is replaced by 1 — 6 X . 

The sum of both terms in Eq. (3.71) produces the matrix element of the scattering matrix on-energy-shell 
which is independent of A. Our complete on-shell result in the effective theory is equal to the well known 
Feynman result for the one boson exchange scattering amplitude. 

There is an important property of the second order calculation above which is worth a separate note. 
When the hamiltonian width in the mass difference becomes small the effective meson emission can no longer 
occur. Thus, the effective theory describes fermions interacting by potential forces. The potentials are given 
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by factors .f (zf 2 ,\).922A in Hti\- The form factors f(z 2 2X ) are known off-energy shell. g 2 2A contains also the 
inner similarity factors which force the intermediate boson to form a sufficiently high invariant mass state 
but if the width is small enough these factors are equal 1. f{z 2 2:K )g22\ is the generalized relativistic potential 
term that equals Yukawa potential in the nonrelativistic limit. Thus, we have accomplished a derivation of 
the boost invariant potential term order g 2 in the effective Yukawa theory. 

The nonrelativistic Yukawa theory is obtained when the width A is such that the allowed energy transfers 
are much smaller than the effective fermion masses. This condition limits only the relative motion of the 
effective fcrmions. It does not limit their total momentum which can still be chosen arbitrarily by taking 
advantage of the boost invariance. The reduction of the fermion dynamics in Yukawa theory with small A 
to the Schrodinger equation in second order perturbation theory is further discussed in Ref. |2(|. Initial 
studies of 4th order similarity in a Yukawa model are given in Ref. |27j] . 

Note that in the energy diagonal part of the effective potential as well as in the on-energy shell scattering 
amplitude the outer similarity factor equals 1 independently of the size of the momentum transfer. In other 
words, one cannot see the outer similarity factor in the physical scattering amplitude order g 2 and the only 
signs of the effective nature of the potential are the form factors in the interaction vertices. 

If we use the interaction gipi^ftpcj) instead of gipiprf' in writing the initial hamiltonian of Eq. (3.3), the 
resulting effective potential corresponds to the one-pion exchange between nucleons. Since the formalism is 
not limited to the nonrelativistic domain of the fermion momenta or to the lowest order perturbation theory, 
one can investigate this type of potentials in a wide range of applications in meson-baryon and quark-pion 
physics. 

There exists a possibility that the similarity flow of hamiltonians may lead to growth of coupling constants 
for small width. The outer similarity factor reduces the strength of the effective interactions when A decreases. 
Effective hamiltonians with small width may have the same bound state eigenvalues as hamiltonians of similar 
structure with large widths and small couplings if the effective coupling constants become large for small 
widths. The range of coupling constants requires investigation in order to establish if the size of coupling 
constants required in the meson-nucleon phenomenology can be explained this way. 



The fermion-anti-fermion interaction order g 2 satisfies a differential equation which is analogous to Eq. 
(3.40) but more terms appear. The operator subscripts must distinguish fermions and anti-fermions and 
one has to include terms which result from the annihilation channel. The fermion-anti-fermion effective 
interaction term is 



Note the change of order of the spin numbering and momentum assignments in comparison to Eq. (3.42) 
for fermions. The new order results from the operator ordering including anti-fermions as defined in Section 
2. a. Momenta fej and fc 2 are used for fermion and k 3 and for anti-fermion operators with even subscripts 
for annihilation operators and odd subscripts for creation operators. 

There are three terms contributing to the derivative of <?iniA with respect to A: one due to the annihilation 
channel and two due to the exchange of a boson. One of the latter two contributions results from the emission 
of the boson by the fermion and absorption by the anti-fermion and the other one from the emission by the 
anti-fermion and absorption by the fermion. In each of the terms there are two similarity functions with 
different arguments. We have 



Fermion-anti-fermion interaction 




(3.73) 



dguiix 



1) 8{x — 

^— + 5 , 3?'32'"3512^33' , 3534 7 ; 



dX 



- |S , 2 r 2 ir 2 5i 2 r 24 r 2 534 
The inner similarity factors are 



(y - x)P 



9{y-x) 




} 



guiu 2 gviv 3 . 



(3.74) 



Si = /(4) 



[-/(4)]' h/(4)]' 



/(4)« 



(3.75) 



AE i2 AE a 
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Equation (3.75) is similar to Eq. (3.44) (the subscript A is skipped for clarity). The second subscript of the 
arguments of the similarity function / denotes the vertex, i.e. 1 stands for the vertex where the boson was 
annihilated and 2 stands for the vertex where the boson was created. In Eq. (3.74), the fermion regularization 
factors first subscript is the same as the corresponding inner similarity factor subscript (i.e. the subscript 
of S) and the second subscript is the same as the corresponding fermion momentum subscript. The boson 
regularization factors are distinguished by the subscript 5 following the convention from Eqs. (2.50). Their 
first subscript is also the same as the corresponding inner similarity factor subscript. Last two subscripts of 
the boson regularization factors equal subscripts of the fermion momenta from the vertex where the boson 
regularization factor originated. Arguments of the regularization factors have the same subscripts as the 
regularization factors themselves, i.e. T{ = r(eei/ A 2 ). The daughter energies in the arguments are calculated 
according to the rules given in Eqs. (2.47) to (2.50). We give the results below for completeness. The same 
arguments will appear in all theories of physical interest. 



en = 



n ±2 + m 2 



ei3 



ei2 



e2i 



C2512 



k ±2 +m 2 
1-x 



p ±2 +m 2 

V 

p ±2 +m 2 
/4J + m 2 

X212 

k ±2 -A- II 2 
_ K 212 ' M 

1 - £212 



X 

X212 = ~ ■ 
V 



K 212 = K± ~ X212P 1 - 



e24 = 



e2534 



,±2 
l 234 



X234 



,±2 
' 2.", J. 



^234 = 



£234 

1-2/ 



^234 



e32 



-P + X 2 34K 



v 312 



C3512 



£312 

K ±2 4- II 2 
K 312 ' M 

1 - £312 

y 

X312 = - ■ 

X 



K 312 = P 1 ' ~ x 312 ■ 
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e 33 = K ™ l + m2 . (3.78.e 

X334 

e3534 = ^4±^. (3.78./ 

1 - £334 

1 — X 

2:334 = r. • (3.78.5 

1-2/ 

ft 334 = -k^ + £334/^ • (3. 78. ft 

Arguments of the similarity functions and energy denominators which appear in Eq. (3.75) are calculated 
according to the rules given in Eqs. (2.12) to (2.19) and (2.24) to (2.26). The results are universal for all 
one-particle-exchange two-particle interactions and are given below for completeness. 

AM\ X = (fci + k 3 )l - (k lm + k 3m ) 2 

= yu 2 — en — ei3 . (3. 79. a 

EA^n = —AMi! + 2fi 2 . (3.79.6 

AE n = AM 2 n /P + . (3.79.C 

AM 2 12 = (k 2m + k 4m ) 2 - (k 2 + k 4 )l 

= ei4 + ei2 - [i 2 . (3.79.d 

Y,M\ 2 = AM\ 2 + 2^ 2 . (3.79.e 

AE 12 = AM 2 12 /P+ . (3.79./ 

AM\ X = (fc 2 534p + fc 4 m) 2 - fcL 

= e 2 534 + e 2 4 - m 2 . (3. 80. a 

T,M 2 21 = AM 2 21 + 2m 2 . (3.80.6 

AE 21 = A7W| 1 /(1 - x)P+ . (3.80.C 

AM 2 22 = k 2 m - (fc 2 5i2 M + fci m ) 2 

= to 2 - e 2 5i2 - e 2 i . (3.80.G? 

Y,M 2 22 = -AM 2 22 + 2^ 2 . (3.80.e 

AE 22 = AM 2 22 /yP+ . (3.80./ 

AM 2 31 = (/c 3 5i2 M + k 2m ) 2 - k\ m 

= e 35 i 2 + e 32 - m 2 . (3. 81. a 

Y,M 2 31 = AM 2 31 + 2m 2 . (3.81.6 
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AE 31 = AM 2 31 /xP+ . (3.81. c) 

AM 2 32 = k\ m - (k 3534f , + fc 3m ) 2 

= m 2 - e 35 34 - e 33 . (3.81. d) 

ZM 2 32 = -AMj 2 + 2fi 2 . (3.81.e) 

AE 32 =AM 2 32 /(l-y)P+. (3.81./) 

In all cases, the arguments of the similarity functions are given by Eq. (2.18), i.e. = AMf / (EMf + A 2 ) 
for all subscripts appearing in Eq. (3.75). 

The same reasoning is used to integrate Eq. (3.74) as in the case of Eq. (3.65) for the fermion-fermion 
interaction. For the similarity function / approaching the ^-function with uo = 1/4 in Eq. (2.20), we have 
f(zf) = 9(X 2 + 2m 2 - |A.M 2 |) with mf = [i 2 in the first, and m 2 = m 2 in the second and third inner 
similarity factors in Eq. (3.74). 

Integration of Eq. (3.74) gives 

SiIiia = 9iiiie + ci guiv 3 gv 4 u 2 r n r 13 r 14 r 12 
+ [C2 r 21 r 2512 r 24 r 2534 9(y - x) + c 3 r 32 r 3512 r 33 r 3534 9(x - y)] gu 4 u 2 gv 4 v 3 , (3. 82. a) 

where the coefficients arc, 



#12-11#12-A . #11-A#11-12 ., , 

Cl = lAMlT + \AM 2 1 ' (3 ' 82 - 5) 



^^22-21^22-A (1 — x)^21-A^21-22 ,„ „„ , 

° 2 - (y-x)\AM 2 22 \ + (y-x)\AM 2 21 \ ' 

_ (1 - y)9 3 2- 3 i0 3 2-\ x9 3 i-\8 3 i^ 32 
C3 - (x-y)\AM 2 2 \ + (x-y)\AM 2 31 \- 

The symbols for ^-functions have the following meaning. 9i-j = 9(\AM 2 \ — \ AM 2 \) and 9i-\ = 9(\AM 2 \ — 
2m 2 — A 2 ) with m 2 equal pi 2 in c\ and m 2 in c 2 and c 3 . 

The next step is the construction of the interaction F\ [Ghha] from 5mu of Eq. (3.73) using Eqs. (2.8) 
and (2.9). 

Then, one has to find out if matrix elements of F\ [Giiiia] between finite free invariant mass states have 
a limit when e — > 0. Stated differently, one checks if the existence of the limit requires the initial value of 
ffiiiie to differ from zero to cancel potential divergences in the limit. Following the same steps as in the case 
of Eqs. (3.66) and (3.70), one can check that no divergences arise. Therefore, gmu = 0- 

The final answer for the effective fermion-anti-fermion interaction is 

#iiiia = F\ [Giiiia] 

= J[P]j5+ J2 /MM5llllA/(2 1 2 niA) & Lp+ KCTl 4(l-^)P- K( T3^(l-2/)P-P^ 6 A 1y P+pa 2 , (3.83.o) 



where 



SiiiiA = ci guiv 3 gv 4 u 2 + [c 2 9(y - x) + c 3 9{x - y)] gu x u 2 gv 4 v 3 . (3.83.6) 
In terms of the fcrmion momenta, 



9(s — 3^i 2 - A 2 ) /ooo , 

ci = -i — ? '- , 3.83.C 

3 - [l l 



with s = max(M 2 3 ,M 24 ), and 
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c 2 % -x) + c 3 6(x -y) = Ua - bUa - x + »±2*±± . (3 .83.d) 

a b 

with a = n 2 -ql 2 , b = ^ 2 -q\^ Q a -b = l-#6- = 0{m xy a-mi- xl - y b), 9 a -\ = [m xy a - \x - y\(2m 2 + A 2 )] , 
8i,-\ = 9 \m\-xi-yb — \x — y\(2m 2 + A 2 )] , m xy = max(x, y) and mi- x i- v = max(l—x, 1—y). The argument 
of the outer similarity factor in Eq. (3. 83. a), i.e. ZiniAj is equal to z 2 2\ from Eq. (3.70.h). Note that Eqs. 
(3.83.a-d) provide the generalization of Eqs. (3.70.a-h) to the case of effective interactions of distinguishable 
fermions. 

When A 2 is reduced below 4m 2 — 3/x 2 , the internal similarity factor in the annihilation term stays equal 
1 independently of the value of A. The effective interaction term provides the full contribution of the 
annihilation channel to the fermion-anti-fermion scattering amplitude of order g 2 . The fermion-anti-fermion- 
boson term in the effective hamiltonian which could contribute acting twice in the scattering is zero for such 
low values of A 2 because the mass gap between the boson and the fermion pair is larger than A allows. 

The internal similarity factor in the exchange term becomes equal 1 independently of A only when A 2 
becomes smaller than —2m 2 + 2m fi + /i 2 . The lower bound on A 2 is — m 2 — (m + ix) 2 (see the discussion of 
Eqs. (2.18) and (2.20)). In the lower bound region, the effective boson emission and absorption vanish and 
the exchange interaction term provides the full scattering amplitude due to the one-boson exchange. The 
amplitude is equal to the standard result on-shell where the outer similarity factor equals 1. 

If the boson mass is much smaller than the fermion mass the fermion energies for small momenta are 
quadratic functions of momentum while the boson energy is a linear function of momentum. Therefore, for 
sufficiently small momenta, the boson energy is large in comparison to the fermion kinetic energies and their 
changes. Thus, the one-boson-exchange interaction is mediated by a relatively high energy intermediate 
state. Consequently, it is contained in a potential term in the effective hamiltonian. 

For small A, the effective hamiltonian contains potentials which are equal to standard scattering ampli- 
tudes in the Born approximation. The potentials differ from the Born amplitudes off-shell in a unique way 
which is dictated by principles of the hamiltonian quantum mechanics and the similarity renormalization 
group: the outer similarity factor reduces the strength of the interaction off-energy-shell. In the light-front 
dynamics the off-shellness is measured in terms of the free invariant mass. 

Outside the lower bound region for A the scattering amplitudes obtain also contributions from the effective 
interactions which change the number of bosons by one in the transition through the intermediate states. 
Analysis of Eq. (3.71) in application to the fermion-anti-fermion scattering follows the same steps as for the 
fermion-fermion scattering in the previous Section. The resulting on-shell scattering amplitude is independent 
of A. The amplitude is equal to the well known perturbative result in Yukawa theory to order g 2 . 



3.c QED 

This Section describes calculations of the effective mass squared term for photons, the effective mass 
squared term for electrons and the effective interaction between electrons and positrons in QED. The calcu- 
lated terms are order e 2 . 

The initial expression which we use to calculate the renormalized hamiltonian of QED is obtained from 
the lagrangian C = — jF^F^ + ^>{i If) — m)4> by the procedure of evaluating the energy-momentum tensor 
and integrating T H over the light-front. |23] We have 



H. 



QED 



dx d 2 x A 



-d ±2 + m 2 
2id+ 



Ipr, 



Iap ±2 a 0v 



e 2 ipn 



4o 



2id+ 



+ TTlpmT'tpii 



(id+y 



(3.84) 



x + =0 



where ip m is a free fermion field with mass m and Aq is a free massless photon field with A J = 0. 

We replace fields ip m (x) and Aq (x) for x + = by the Fourier superpositions of creation and annihilation 
operators, we order the operators in all terms and we drop terms containing divergent integrals which result 
from contractions. This is done in the same way as in the Yukawa theory but more terms need to be 
considered. Then, we introduce the regularization factors. 

The ultraviolet regularization factors are already familiar and the same as in the Yukawa theory. The 
additional regularization is required due to the infrared singularities. Photons have diverging polarization 
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vectors when their +-momentum approaches 0. The corresponding seagull term, i.e. the 5th term in Eq. 
(3.84), is diverging too. We introduce the infrared regularization factors (1 + 6/x)~ 1 as described in Section 
2.b. 

We also introduce a photon mass term p% — /i 2 by adding it to — d ±2 in the second term in Eq. (3.84). 
A fixed value of fj,$ leads to the conclusion in perturbation theory that the photon eigenstates have masses 
equal to fig when the charge e approaches 0. Therefore, we will be forced to consider the limit fi$ — > in 
order to discuss physical photons to order e 2 . Also, the nonzero mass squared term for photons leads to 
additional divergences when S — > and the limit of /i$ — > removes those. 

The infrared finiteness of QED suggests that physical results in our approach should be independent of 
fi$ when it is sufficiently small. We introduce the photon mass /is and investigate the limit fig — * 0. The 
second order calculations in this paper lead to results which are independent of ji^ when it tends to zero. 



Photon mass squared 

The same procedure from Section 2 which led to Eqs. (3.4) and (3.5) in Yukawa theory leads in QED to 



Gi photon x = ^2 J W ~ k + r ' X a lo- a ko- ■ (3.85) 

No correction arises to the term k 1 - 2 /k + because our regularization preserves the kinematical symmetries of 
light-front dynamics. 

A new feature in comparison to the Yukawa theory is the polarization of photons. With the kinematical 
symmetries explicitly preserved, only terms diagonal in the photon polarization emerge. For example, terms 
proportional to k l e l ai k J e J a2 with o\ ^ <y% cannot appear because the regularization and similarity factors 
do not introduce dependence on the photon momentum. Note that such terms are allowed by the power 
counting. || 

The net result of the photon self-interaction is an effective photon mass squared term which is independent 
of the photon momentum but varies with the effective hamiltonian width A. One obtains more complicated 
results for the effective photon free energy if the regularization or similarity factors violate kinematical 
symmetries of light- front dynamics. || 

The dependence of /u? on A is determined to order e 2 by the equation 



~dX S =e J [XK] —A M^ 2 r ^ K) (3 - 86) 

and the initial condition at A = oo. However, it is also sufficient to know the effective mass squared at any 
single value of A to determine its value at other values of A using Eq. (3.86). The initial condition at A = oo 
is distinguished only because it provides connection with standard approaches based on the local lagrangian 
for electrodynamics. 

In Eq. (3.86), M 2 = (k ±2 + m 2 )/x{\ - a:), AM 2 = M 2 - [i\ and Y.M 2 = M 2 + nj so that z x = 
(M 2 -fi 2 s )/{M 2 + fi 2 s + \ 2 ). In the limit of Eq. (3.8), we have / 2 (z 2 ) = 6{\ 2 + 3$ -M 2 ). In fact, Eq. (3.86) 
is free from infrared singularities and we could skip the introduction of /i 2 by letting it go to zero at this 
point. However, the systematic approach defines the hamiltonian of QED including the infrared regulator 
mass for photons and we can keep it here for illustration. The regularization factor r e (x,K) in Eq. (3.86) 
is the same as in Eq. (3.9) in Yukawa theory because only fermion regularization factors enter Eq. (3.86), 
according to Eqs. (2.47) to (2.49), and these factors are the same in both theories. 

Evaluation of the spin factor gives 



e 2 / M ^_A_AZ 1_ _ rt[X} K)i ( 3. 87) 



M_ 2 [t df(zl) 2M 2 -4k ± 2 
dX - e J [xk\ dX M2 _ ^ 

which is the QED analog of Eq. (3.6) from Yukawa theory. 

Integration of Eq. (3.87) is carried out through the same steps as in the Yukawa theory. We can use Eq. 
(3.22) to calculate the effective photon mass squared fi\ knowing its value /x 2 , at some value of A = Ao. 

The value of /i 2 , is found by requesting that the effective hamiltonian eigenvalues for photon states contain 
the physical photon mass fi, expected to be 0. However, solving the eigenvalue equation to second order in 
the coupling constant e through the same steps as in the case of mesons in Yukawa theory in Eqs. (3.23) 
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to (3.25), leads to the physical photon mass fl = n$. /is is small but finite and it is considered, in terms of 
powers of e, to be of order e° = 1 when e — > 0. 

In the ^-function limit for the similarity function / with uq = 1/4 one obtains 

" s = "* + £ jf ix f + 3 "« - M "> + * 4 > • (M8) 

Thus, at other values of A, the effective photon mass squared is 

q pi poo 

47t7 dX J 



This result naturally depends on the infrared regularization parameter /j, 2 but no singularity appears when 
this parameter is set equal to zero. For A 2 + 3/x 2 < 4to 2 , where 4m 2 is the lowest possible free invariant 
mass squared for the two intermediate fermions, the photon mass is independent of the hamiltonian width 
A and it equals p 2 . For larger values of A 2 , the effective photon mass grows with the width A so that its 
larger value compensates effects of the interactions which become active for the larger width. The net result 
is that the photon eigenstates have eigenvalues with masses squared equal [i 2 independently of A. Finally, 
the result favored by experimental data is obtained in the limit \x\ — > at the end of the calculation. 



Electron mass squared 

Electron and positron self-interactions through emission and reabsorbtion of transverse photons result 
in the fermion free energy terms of the form exactly the same in QED as in Eq. (3.27) in Yukawa theory. 
However, the effective mass of electrons and positrons depends on the width differently than in the case of 
Yukawa theory. Instead of Eq. (3.28), one obtains now 



XK 



d f 2 ( z \) Um ° k 4<r(An + TO) f^-Umak 



M 2 - to 2 



■r eS (x,K), 



(3.90) 



where k = (k -p)o, P + = xk + , p 1 - = xk 1 - + Kr 1 , M 2 = (to 2 + k 2 )/x + (p 2 + k 2 )/(1 - x), AM 2 = M 2 - to 2 , 
T,M 2 = M 2 + m 2 and z\ = AM 2 / (T^M 2 + A 2 ). The regularization factor of Eq. (2.49) for the intermediate 
particles and the infrared regulator for the intermediate photon, as given at the end of Section 2.b, imply 



r eS (x, k) 



-M 



e \ 2 k 2 + TO Z K A + fig 

x 1 — X 



n 2 „2 



1 -2 



f 



1-X 



A 2 "" 1 ' (a 2 ) 

The sum over photon polarizations in Eq. (3.90) produces the well known expression 

k a g +fj + g +a k & 



' / fecr ka a 



k+ 



(3.91) 



(3.92) 



and the spin factor in Eq. (3.90) is 



777,(7 A; 



_ a0 , ~k a g +ti + g +a ¥ 
9 ' " k+ 



(1 - x) 2 m 2 + k 



1 + x 2 



(3.93) 



The new feature of this expression, in comparison to the Yukawa theory, is the divergence for x — > 1, i.e. 
where the photon longitudinal momentum approaches 0. 

The rate of change of the electron mass term versus the effective hamiltonian width in the ^-function 
limit for the similarity function with uq = 1/4 is 



d\ 2 wit j 
where u — k 2 /x(1 — x), 



,1 ,oo m 2 2{l-xY 

= JL dx du <5(3to 2 + A 2 - M 2 ) x —„ j^^r eS (x, K ), 

4tt J J M 2 - Tfi 



(3.94.a) 
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and 



M 2 = u + — + . (3.94.c) 

x 1 — x 

The divergence structure of the effective electron mass in Eq. (3. 94. a) is obscured by the fact that the 
whole effective mass term is merely a number dependent on A while three cutoff and regularization parameters 
appear in the integral: e, S and [ig. The only available condition is that the effective electron mass should 
have a limit when e — > 0. However, this condition has to be satisfied without generating divergences in the 
physical electron mass (i.e. in the electron eigenvalue energy) when we remove the infrared regularization. 
Since other contributions to the physical electron mass may diverge as 6 or [i$ tend to 0, and only the sum 
is finite in the limit, one needs to keep track of the infrared structure in defining the e-independent (i.e. 
ultraviolet finite) part of the counterterm. 

The divergences due to e — > 0, 5^0 and fj,g ~ * 0, are not resolved in the single mass constant. Many 
elements of a complete analysis overlap in producing the final answer and many simplifications are possible. 
We will proceed in this Section with a simplified analysis. A more extended analysis will be required for 
other hamiltonian terms where the outcome of the procedure is not reduced to finding only one number in 
the effective interaction. For example, in the electron-positron interaction term the external momenta of 
fermions introduce a whole range of additional parameters. That case will be illustrated in the next Section. 

The ^-function under the integral on the right-hand side of Eq. (3. 94. a) forces Ai 2 to be equal A 2 + 3m 2 . 
The smallest possible value of Ai is m + [is and the negative value of A 2 = (m + fig) 2 — 3m 2 = — m 2 — (to + 
Us) 2 + 2[(m + us) 2 — to 2 ] is required to reach this lower bound (see the comments about Eq. (2.20)). Below 
this bound the right-hand side of Eq. (3. 94. a) vanishes, no emission or absorption of photons by electrons 
is possible and the effective electron mass stays constant. The smallest possible value of A 2 allowed by Eq. 
(2.18) is — m 2 — (to + lis) 2 - The difference between these bounds vanishes when the photon mass goes to 
zero. 

In the next Section about electron-positron interaction we will also consider the case of the infinitesimal 
u (see the discussion below Eq. (2.20)), which leads to 6(X 2 + m 2 — M 2 ) in Eq. (3. 94. a), instead of 
<5(3to 2 + A 2 - M 2 ). 

For A 2 > (m + fig) 2 — 3m 2 , the right-hand side of Eq. (3. 94. a) is positive. Therefore, the effective electron 
mass squared term grows together with the width of the hamiltonian. This growth combines with the growing 
negative contributions of the corresponding effective transverse photon emission and reabsorption so that 
the physical electron mass is independent of the hamiltonian width. 

For A 2 close to the lower bound, u must be close to and x <~ 1 — [Jbg/{m + fig). But M 2 is limited 
and determined by the value of A. Quite generally, as long as Ai remains limited, u and fi 2 /(l — x) are 
limited. The invariant mass denominator in Eq. (3. 94. a) equals A 2 + 2to 2 . The denominator is small only 
when A 2 is close to the lower bound. Then, the integration range is small too. The denominator, when 
expressed in terms of A, can be pulled out and put in front of the integral. The integration over u sets 
u = u(x) = A 2 + 3m 2 — m 2 /x — /i 2 /(l — x) and forces the condition u(x) > 0. This condition implies the 
following limits on the integration over x, provided A 2 + 3m 2 > (m + [is) 2 since otherwise the integral is 0. 

£0 — Ax < x < x + Ax , (3. 95. a) 

-=K i+ ot)' <3 ' 95 - 6) 

Ax = ^x 2 - . (3.95.C) 

Within these limits, u(x) varies from the minimal value of at the lower bound xo — Ax through a maximum 
of A 2 + 3to 2 — (to + fis) 2 at x = m/(m + [is) to the minimal value of again at the upper bound xq + Ax. 
In the case of infinitesimal u , one replaces A 2 + 2to 2 in these formulae by A 2 . Eq. (3. 94. a) reads 
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where 



ad{3m 2 + A 2 - (m + ^ s ) 2 } 
47r(A 2 + 2m 2 ) 



xo — Ax 



,2(1 



u(x) 



2(1 -+a; 2 )1 



r e s(x, u(x)) , (3. 96. a) 



r e5 (x,u(a;)) 
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1-x 



1 



(3.96.6) 



The upper limit of integration over a; for /z 2 much smaller than m 2 and 2m 2 + A 2 , equals 1 — /i 2 / (2m 2 + A 2 ) 
and approaches 1 when ^5 ^ 0. For x close to 1, the factor (1 — x)^ 1 in the square bracket of the integrand 
is large and leads to a logarithmic dependence of the integral on the upper integration limit. The logarithm 
would become infinite for [15 — > if 6 were equal 0. Therefore, the limit /ig — > is sensitive to the presence 
of the regularization factor with S ^ 0. For the finite /i«5, the region of x — > 1 is regulated by e and can be 
considered an ultraviolet limit. A counterterm to the diverging e dependence could remove the divergence 
due to x — > 1. Then, a separate cutoff parameter 8 would not be needed but the resulting terms would 
diverge for fj, 2 — > 0. For finite /i 2 , the derivative of the electron mass with respect to A is finite. For fi 2 = 0, 
the integrand in the region x ~ 1 is regulated solely by the infrared regularization factor [1 + 5/(1 — x)]~ 2 
since the upper limit of integration over x is equal 1. 

For finite A 2 all three terms in M 2 , i.e. u(x), m 2 /x and ^ 2 /(l — x) are limited. Therefore, for finite [is, 
one can take the limit e — > in the integrand. The factor [1 — 6/(1 — x)]~ 2 remains and additionally cuts 
off the integration region at x ~ 1 — 5. The ratio of 6 to n 2 /(2m 2 + A 2 ) determines the size of contributions 
obtained from the upper range of integration over x. For finite fixed values of /is, one can take the limit 
6^0 and log [(2m 2 + \ 2 )/fi 2 ] appears in the answer. 

The right-hand side of Eq. (3. 96. a) contains terms which behave for large A as a constant, as A~ 2 with 
factors of logarithms of A and as functions vanishing faster than A~ 2 . One integrates Eq. (3. 96. a) over A 
from Ao to infinity in order to express the effective electron mass squared term at A , denoted m 2 ,, in terms of 
the initial m 2 . Clearly, the integration over A would diverge without the regularization factor which depends 
on e. The integration produces terms behaving as e _1 , loge and terms convergent in the limit e — > 0. m 2 in 
the initial hamiltonian must be supplied with a counterterm to subtract the diverging e-dependent terms in 
the effective hamiltonians. 

In summary, the infrared divergence due to us — * and 6 — > appears in the derivative of m 2 with 
respect to A. Therefore, even if one requests that the electron mass term is finite at some value of A, the 
effective masses of electrons in the neighboring hamiltonians with even slightly different widths will diverge 
when us and 6 approach 0. We have to abandon the requirement that the effective electron mass term at 
any value of A remains finite when the infrared regularization is removed. The effective masses diverge in 
the limit fj,$ — » and 6 — ► 0. The only condition we can fulfill through the ultraviolet renormalization is 
that the effective electron masses for finite A are independent of e. 

Mathematical details of the effective electron mass term calculation are more complicated than in the 
Yukawa theory because the infrared regularization parameters are present. Otherwise, the calculation is 
essentially the same and we skip the description here. We only stress that the counterterm and the effective 
masses of electrons and positrons depend on the infrared cutoffs and they diverge when the cutoffs are being 
removed. 

Thus, the effective electron mass squared term in the limit e — > is 



m 2 2(l-x) 2 »=» 2(l+x 2 ) 2 

m 2 =m 2 + e 2 J [xn] [f 2 (z 2 ) f(z 2 o )} ~ (l + + o(e±) . (3.97) 

The finite term m 2 , has a limit when e — > 0. Its dependence on the infrared regularization is not displayed. 
m,Q is found from a suitable renormalization condition. 

The natural condition to be satisfied by m 2 , is that the effective hamiltonian of some width A has the 
electron eigenstates with eigenvalues equal (p ±2 + m 2 )/p + , where p denotes the electron momentum and m 
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is the physical electron mass. The eigenvalue equation for electrons can be solved in perturbation theory in 
the same way as for bosons in Yukawa theory in Eqs. (3.23) - (3.26) and fermions in Eqs. (3.35) - (3.39), or 
in QED for photons in Eqs. (3.88) - (3.89). One obtains the condition 



m 2 



2 2(l-x) 2 . J 2(l+s 2 ) 2 

= ml-e>J[ XK ]f( Z l) * M2 *~? X ~ X ( 1+ i~j + (3.98) 

and one can calculate mo from this condition. Consequently, 

22(1^ f 2(l+x 2 ) 2 

ml=f^ + e' j[xn]f{zl) ' ^ 2 ^ (l + + o( e 4 ) , (3.99) 

and m 2 = to 2 + o(e 2 ). 

The physical electron mass is independent of the infrared regularization because the regularization de- 
pendent m\ and the effective emission and absorption of photons combine to the regularization independent 
result. 



Electron-positron interaction 

Calculation of the effective electron-positron interaction to order e 2 is of interest as a way to derive the 
Coulomb force in quantum electrodynamics - this interaction is responsible for the formation of positronium. 
Also, effective interactions between quarks and anti-quarks in QCD have a number of similar features and 
the QED calculation provides an introduction to the QCD case. 

Generally speaking, the QED calculation of the effective electron-positron interaction proceeds in the 
same way as in the case of fermion-anti-fermion interaction in Yukawa theory except for three new elements. 

The first is that photons have polarization vectors which enter in the vertex factors and introduce addi- 
tional dependence on the exchanged photon momentum. This dependence leads to infrared divergences for 
small longitudinal momenta of exchanged photons. 

The second feature is that the infrared divergences require additional regularization factors. We use the 
parameter <5 and the photon mass squared us ^ 0. The limits S — ► and ^5 ^ are generally understood 
as to be taken at the end of a calculation of observables and not in the effective hamiltonian itself. However, 
it may also be possible to take the limits in matrix elements of the hamiltonian between states which do not 
induce infrared divergences, i.e. do not involve small x photons. 

The third feature is that the one-photon exchange interaction needs to be combined with the 5th term 
from Eq. (3.84) to obtain the standard results for the electron-positron scattering in the Born approximation. 
The 5th term from Eq. (3.84) provides the initial condition for the renormalization group flow of the effective 
hamiltonians. In order e 2 this term is only supplied with the outer similarity factor by the operation F\. It 
does not change in the flow beyond this factor because it is order e 2 itself. The initial condition provides a 
contribution which is needed to obtain the Coulomb potential. This is a different situation than in Yukawa 
theory where no four-fermion seagull interactions appeared and the one-meson exchange interaction was 
sufficient to produce the Yukawa potential in the effective hamiltonians of small widths. 

The effective electron-positron interaction term has the analogous structure as the fermion-anti-fermion 
interaction in Eq. (3.73), i.e. 

Sum = /[ p ]^+ Yl f[ XK ]lyp}9mi\btp +Kai dl 1 _ x)p _ Ka3 d {1 _ y)P _ pa4 b yP+pa2 . (3.100) 

(71(72(7^(74 

The coefficient function guiiA °f order e 2 satisfies the differential equation 

"0-5 
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-£3^32^3512^33^3534^1/5^/5— ——— 2jeui j4 5<Ts u 2 ev4 fk 5 a 5 v 3 , (3.101) 

(x y)f ^ 

which is the QED analog of Eq. (3.74) from Yukawa theory. Notation is the same as in Eqs. (3.74) to (3.81) 
with the exception that fi 2 is replaced by fj, 2 . The new elements are the infrared regularization factors of 
Section 2.b, i.e. r 3 / 5 = r(k^5/k^), r 2 / 5 = r(k 2 S/k£), r\/§ = r(k 1 S/k^'), r 4 / 5 = r{k 4 5/k§), and the photon 
polarization vectors. The sum over photon polarizations gives 

Y, £ t°AL - -9 af3 + 509 t ■ (3 - 102) 

The terms proportional to the four-vector fcso can be rewritten using the Dirac equation for free fermions of 
mass to. For example, in the second term on the right-hand side of Eq. (3.101) we have 
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Using similar relations for all vertex factors involved one obtains 
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(3.103) 



(3.104) 



We use the notation = (fcj + kj) 2 and q 2 j = (ki — kj) 2 . 

Integration of Eq. (3.104) proceeds in the same way as in the case of Eq. (3.74). The initial condition 
at A = oo includes the seagull term. 
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The result of the integration of Eq. (104) is 
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The coefficients c\, c 2 and c 3 are given by universal Eqs. (3.82.b) to (3.82.d) which were derived already in 
Yukawa theory, with the replacement of the meson mass by the photon mass, /j, 2 = (/%. 

This result will be now analysed term by term for illustration in the electron-positron scattering in 
second order perturbation theory. The first term in Eq. (3.71) for the T-matrix calculated to order e 2 has 
matrix elements equal to the matrix elements of the effective interaction from Eq. (3.106). When evaluating 
the ^-matrix elements, one considers configurations where the free energy of incoming fermions equals the 
free energy of the outgoing fermions. This configuration selects the energy-diagonal part of the effective 
interaction: s 13 = ,s 24 = s and q 2 2 = q 2 4 = q 2 . Also, the external similarity factor which appears in the 
effective hamiltonian as an additional factor to gniiA equals 1. Thus, gnn\ in the energy diagonal part can 
be viewed as the scattering amplitude. It simplifies in the energy diagonal part to 

Sum = 9 miT rterm ~ cie«i7 M «3 ev^u 2 

+ [-cis + 1] — 

- [c 2 r 2 5i 2 r 2 5 34 r 3/5 r 2/5 (% - x) + c 3 r 35 i 2 r ? , 53i r l/b r A/b 9(x - y)] eu 1 j fi u 2 ev 4 1f,v 3 
- { [c 2 <? 2 + l] r 2 5i 2 r 2 534r 3/5 r 2/ 5% - x) 
,r 2 , ,1 ai u eui7+u 2 ew 4 7 + W3 

+ |c 3 <7 + 1J r 35 i2r 353i r 1/5 r i/5 0(x -y)) — ( x _ y yp+2 — ■ (3.107) 

We have removed the regularization factors which equal 1 for fermions with finite free energy. 

In the limit (is — ► and for the cutoff A close to —2m 2 , a number of simplifications occur. We display 
the result for the case where the intermediate photon momentum fraction \x — y\ >> (x, y, 1 — x, 1 — y)6, 
i.e. when the photon momentum is not negligible in comparison to the fermion momenta. In this case, the 
infrared regularization factors for the photon which are still kept in Eq. (3.107) equal 1 and 

_ n counterterm "17^3 ^47^2 

Slim - flillle 

+9 [\q 2 \ -\x- y\(2m 2 + X 2 )/max(x, y,l-x,l- y)} e ^ u ^W* 

-9 [\x - y\(2m 2 + \ 2 )/max(x,y, 1 - x, 1 - y) \q 2 \] ^g^^ffi . (3.108) 

The first term is the potentially necessary counterterm which we have not yet determined. Since the 
remaining terms are not sensitive to e, the counterterm matrix element is equal zero. 

The second term is the well known expression for the electron-positron annihilation channel scattering 
amplitude in the Born approximation. No limits on the fermion momenta appear because the invariant mass 
squared of two fermions is larger than 4m 2 which is the minimal invariant mass difference in the transition 
between the two fermions and a one massless photon state. Since we assume the width A 2 to be small the 
effective hamiltonian contains the full amplitude for transition through the intermediate photon state. 

The third term is equal to the standard second order expression for the electron-positron scattering 
amplitude via one-photon exchange, except for the ^-function factor which forces the momentum transfer to 
be sufficiently large. 

The meaning of this restriction is visible in the T-matrix. The third term contributes to the e + e~- 
scattering amplitude through the first term in Eq. (3.71). The same contribution would originate from the 
second term in Eq. (3.71) if we were using the initial hamiltonian to calculate the T-matrix. In contrast, 
the effective hamiltonian with the small width A limits the effective photon emissions and absorptions to 
small momentum transfers and, therefore, it is not able to provide this contribution through the second 
term in Eq. (3.71). This contribution is then contained in the effective hamiltonian above and comes in the 
scattering matrix through the first term in Eq. (3.71). 

The fourth term is unusual in the sense that it should not appear in the electron-positron scattering at 
all. The fourth term distinguishes the z-axis in its structure and diverges when x — > y. The ^-function factor 
in the fourth term is equal 1 where the ^-function factor of the third term is equal 0. And vice versa, the 
fourth term 9- function equals where the third term 9- function equals 1. 
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The need for the fourth term becomes clear when one recalls that the second term in Eq. (3.71) also 
contributes to the electron-positron scattering amplitude. The relevant contribution comes through the 
effective one-photon exchange which results from the double action of Hi\. Hi\ is given by the operation 
F\ applied to the third term of the QED hamiltonian from Eq. (3.84), (see Eq. (2.11)). The operation 
Fx multiplies the photon emission and absorption vertices by the factor f\. This factor was set equal to 
a ^-function in the current example. Each interaction provides one factor of the 0-function. The resulting 
factor in the second term of Eq. (3.71) is the same in QED as in Yukawa theory in Eq. (3.72), except for 
the antisymmetrization effect which leads to the ^-function which stands in front of the fourth term in Eq. 
(3.108). 

Now, the second term in Eq. (3.71) contains spin factor which is the same as in the last term in Eq. 
(3.104). The g^ v part complements the third term in Eq. (3.108) and produces the full well known one- 
photon exchange scattering amplitude which is free from the ^-function factor. The remaining part provides 
the term which cancels the odd fourth term in Eq. (3.108). Thus, the effective hamiltonian calculated to 
second order in powers of the charge e contains an odd term and 6- functions which are required to compensate 
for the odd contributions to the scattering amplitude from the effective hamiltonian order e acting twice. 
The above analysis of the energy diagonal part of the second order effective hamiltonian explains the role of 
different terms in Eqs. (3.106) to (3.108). 

The analysis also suggests that apparently infrared diverging terms in the effective hamiltonian may 
mutually compensate their diverging contributions in the scattering amplitude on energy-shell. In the current 
example, we see the interplay between the second-order seagull term and the double action of the first order 
emission and absorption of photons. The first order hamiltonian matrix elements diverge when the photon 
longitudinal momentum approaches zero. The second order seagull term compensates this divergence in the 
on-energy-shell T-matrix elements. 

The remaining point to make here is that the result of Eq. (3.106) with the counterterm equal leads 
to the effective light-front hamiltonian version of the Coulomb force in the limit of small A 2 + 2m 2 <C am 2 . 
The key elements in deriving this conclusion are the outer similarity factors and the smallness of a. The 
outline of the derivation is following (cf. Refs. jll| and J26|). 

If only the small energy transfers are allowed by the outer similarity factor, i.e. transfers much smaller 
than the electron mass, then the wave functions of the lowest mass eigenstates of the effective hamiltonian 
are strongly peaked at small relative electron momenta and they fall off very rapidly as functions of the 
relative momentum. This is not true without the outer similarity factor because the function <7iniA alone 
is too large at the large energy transfers and it would produce singular contributions in the large relative 
momentum region making the eigenvalue problem sensitive to the ultraviolet regularization cutoffs. 

Below the width scale the wave functions fall off as dictated by the eigenvalue equation with small a. 
Above the width scale the fall off is very fast due to the similarity factor which justifies restriction to momenta 
much smaller than m and the nonrelativistic approximation for all factors in Eq. (3.106) becomes accurate. 

In the nonrelativistic approximation, q\ 2 = q 2 4 = q 2 and Eq. (3.107) applies. Further, the ^-functions 
in Eq. (3.108) become effectively equal 1 and 0, respectively. The last term is not leading to important 
contribution despite its divergent longitudinal structure because it is canceled by the effective massless 
photon exchange as described earlier in this Section. 

The dominant contributions are provided by the second and the third terms from Eq. (3.108) which are 
well known to have the right nonrelativistic structure for predicting positronium properties in the Schrodinger 
equation. Now, the outer similarity factor becomes irrelevant to the spectrum in the leading approximation 
because the coupling constant is very small (cf. Jll[] and [^6|). In the dominant region the electron velocity 
is order a, the nonrelativistic approximation to the full dynamics produces wave functions with relative 
momenta order am and the outer similarity factor in the effective interaction can be replaced by 1. 

We can use the infinitesimal uo in Eq. (2.20) and replace A 2 + 2to 2 by A 2 (see the discussion below 
Eq. (2.20)). When A is order am and x — y is order a the momentum transfer q 2 is typically order a 2 m 2 
which is much larger than 2(x — y)X 2 in Eq. (3.108). Thus, the ^-function is equal 1 and the second term in 
Eq. (3.108) becomes equal to the standard Coulomb interaction with the well known Breit-Fermi structure 
of the spin factors. This step completes the derivation of the Coulomb potential. The derivation explains 
the effective nature of the Schrodinger equation with the Coulomb potential in the light-front hamiltonian 
formulation of QED. 
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3.d QCD 



The main conceptual complication in hamiltonian calculations in QCD is confinement which is not yet 
fully understood. To order g 2 the manifestation of the confinement problem is the lack of a well defined initial 
condition for the renormalization group flow of the effective hamiltonians. It will require an extended research 
effort to find the class of acceptable initial conditions. For example, the on-mass-shell renormalization 
conditions for the quark and gluon mass terms are questionable. 

We stress the urgent need for the higher order calculations by describing some details of the second order 
calculation of the qq effective interaction. The calculation is similar to the one in QED above, except for 
the option for a different treatment of the last term in Eq. (3.108). Perry [[l2| suggested that the long 
distance part of this term may remain uncanceled in the effective QCD dynamics because of the gluon non- 
abelian gauge interactions. If the last term in Eq. (3.108) would remained uncanceled it could be claimed 
to generate confinement in the light-front hamiltonian approach to QCD. j^] We describe the structure of 
this term in the present approach since it is different than in Ref. fl2]| . The differences result from the 
different definitions of the similarity transformation, boost invariance and not invoking coupling coherence. 
The coupling coherence arguments are replaced by a plain perturbative renormalization condition for quark 
mass terms. 



Quark and gluon mass terms 

Results one obtains from Eq. (2.38) in QCD in second order in G2X can be illustrated by two equations 
for the effective masses of quarks and gluons. The range of widths and coupling constants for which these 
equations can pertain to physics are not known yet. 

We do not write regularization factors in detail. One can write them easily using results of Section 2 and 
previous examples in Section 3. In the case of effective masses a number of simplifications occur as explained 
in the previous Sections. Let us consider an infinitesimal uq in Eq. (2.20) and simplify notation by replacing 
A by A itself. Then, we can write 



dX 



Six = 



df 2 (z x /X 2 )/dX n 
y\2\ — 7; ;=; y2i\ 
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E\\ is the eigenvalue of Q\\ corresponding to the subscript 11. A set of arguments z\ is needed. Namely, 
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m\ and fj,\ axe the effective quark and gluon masses, respectively. Then, 
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(3.114) 



The gluon couples to the quark-anti-quark pairs and pairs of gluons while the quark couples only to the 
quark-gluon pairs. The number of colors above is equal to 3 and the number of flavors to 1. 

These equations are not further studied here for two major reasons. The first one is that we do not know 
the initial conditions to use for such study. The second is that the equations involve two running couplings 
which are not known yet. The third and fourth order calculations are required to find them. 
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The importance of effective mass issue for quarks and gluons is illustrated below by the calculation of 
the small energy transfer effective forces between quarks and anti-quarks. 



Quark-anti-quark effective interaction 

The second order results are similar to QED. For heavy quarkonia one can directly look at Eq. (3.108). 
The only change required is the color SU(3) matrices sandwiched between color vectors of quarks and 
summed over colors of the exchanged gluons. The counterterm is 0. The second term gives the annihilation 
channel potential but the color matrix is traceless and this excludes the single gluon annihilation channel 
from the dynamics of color singlet QQ states. The third term leads to the color Coulomb potential with the 
well known Brcit-Fcrmi spin factors. 

About the last term in Eq. (3.108) it was suggested by Perry Jl2| that a term of this kind may become 
a seed for confining interactions if it is not canceled. A cancelation occurs in perturbation theory when one 
evaluates the model interaction in the QQ sector in the explicit expansion in powers of g to second order 
using the operation R and Eq. (1.3) in the nonrelativistic approximation for the quark relative momentum, 
exactly the same way as for electrons in QED with massless photons. 

However, there are reasons for effective gluons to acquire a gap through the non-abelian terms which are 
absent in effective QED. One can assume that the gluon energies in the QQg sector may be lifted up so that 
the gluons cannot contribute to the model QQ interaction in the way the abelian massless photons can in 
the model electron-positron interaction in QED. As a result of this assumption one obtains the last term in 
Eq. (3.108) acting in the effective QQ sector. 

In the nonrelativistic approximation, the Coulomb term and the term in question are (see Eq. (3.108)) 



g 2 4m 2 



9 



2 



t -\x-y\ 2\ 2 ^f- - \x - y\ 2\* - ? r -2_ ? . (3.115) 

The regularization factors are the same in both terms and they are not displayed. The only effect of 
their presence which matters here is that \x — y\ is limited from below by about 8/2. The initial infrared 
regularization parameter 8 comes in through the initial condition in the renormalization group flow of the 
seagull term. The flow is limited in the second order calculation to the dependence of the outer similarity 
factor on A. The factor 1/2 results from \x — y\/x quar k being limited from below by 8 and the quarks having 
Xquark ~ 1 /2. In fact, the lower bound on \x — y\ is given by 8max(x, 1 — x, y, 1 — y). This is different 
from Ref. fll2f where instead of the ratios of the +- momentum fractions a separate frame dependent scale 
for +-momentum is introduced. 

The third component of the exchanged gluon momentum is q% = (x — y)2m. Thus, we see that the 
uncanceled singular term is represented by the potential which is analogous to the Coulomb potential except 
for that the factor — 1/q 2 is replaced by — 1/q 2 and both terms have mutually excluding and complementary 
supports in the momentum transfer space. 

The seagull ^-function can be rewritten as 9[oj 2 — (\q^\ — uj) 2 — q ], where uj = A 2 /2m. The support of 
this function is two spheres of radius uj centered at q 1 - = and q% = ±uj. The spheres touch each other at 
the point q = q 3 = 0. In this point, g| in the denominator produces a singularity. 

Let us initially consider both terms in Eq. (3.115) as the actual interaction in the model QQ sector, 
i.e. as if they were not affected by the operation R in Eq. (1.3). The Coulomb term works outside the two 
spheres in the cf-space and the singular seagull term works inside. 

In the region of the singularity, both q and q% are small in comparison tow. In this rough analysis one 
can neglect the outer similarity factor 9{\ — |fc 2 — k' 2 \) since it is equal 1 when q = k — k' approaches 0. k is 
the relative momentum of the created QQ pair and k' is the relative momentum of the annihilated QQ pair. 
uj = (A/m)A/2 C A/2 and the spheres have the radius about A/m times smaller than the outer similarity 
factor width in the quark relative momenta, i.e. the relative size of the spheres in comparison to the outer 
similarity factor support approaches when A/m — ► 0. 

Since q 1 - 2 is order q$ in the singular region the divergence when q% — > is logarithmic. The lower limit 
of integration over \q^\ for a given x is given by 2m 8 max{x , 1 — x). However, we assume x = 1/2 + o(g 2 ) 
and we neglect terms of higher order than g 2 in the model QQ hamiltonian. 

The potential resulting from the uncanceled seagull term is given by the following expression (cf. fl2]|), 
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(2tt)3 — qi 

The sign ~ means that the diverging dependence on d is subtracted and the same coefficient stands in front 
of the integral as in the Coulomb potential term. The argument for the infrared subtraction goes as follows. 

If the gluons cannot cancel the last term in Eq. (3.108) they presumably cannot contribute to the model 
quark self-energies either for the same reason. Because the size of the quark mass in the effective hamiltonian 
is unknown one may propose that its value is chosen in the second order calculation in the same way as for 
nucleons in the Yukawa theory in Eq. (3.37) or electrons in QED in Eq. (3.99). A would-be quark eigenstate 
has a finite constituent quark mass when the gluons are allowed to contribute in the whole range of momenta 
from zero up in the eigenvalue equation. This setting is equivalent to the solution Perry proposed for his 
coupling coherence condition for the quark self-energies. Q The argument also illustrates the urgency of 
questions concerning the initial conditions and higher order analysis in the similarity renormalization group 
flow. 

There is nothing wrong with the mass adjustment despite the infrared divergence. We have noticed in 
the previous Section that the arbitrary finite parts of the ultraviolet counterterms can be infrared divergent. 
This time, however, the positive and infrared logarithmically divergent part of the effective quark mass 
term in the model eigenvalue equation for heavy quarkonia remains uncanceled when the transverse gluons 
with 2cj|g 3 | — q 2 > are declared to be absent from the model dynamics. The uncanceled part of the 
effective quark mass term stands in the eigenvalue problem. The point is it can now cancel the diverging S 
dependence in the seagull term which is not canceled because the gluon exchange below A is missing. The 
new cancelation between the incomplete masses and the seagull occurs in the colorless states. It is analysed 
here in the nonrelativistic limit. 

Wc describe the cancelation mechanism in the case of equal masses of quarks. The mechanism is similar 
but not identical to that in Ref . |Q . The infrared divergent mass squared term comes into the quarkonium 
eigenvalue equation divided by x{l — x). But in the second order analysis the mass divergence appears 
only as a logarithmically divergent constant and the ^-dependence is of higher order. The same diverging 
constant with the opposite sign is generated by the seagull term. 

The infrared divergent terms and their cancelation are not directly related to the ultraviolet renormal- 
ization procedure. They appear in the ultraviolet-finite effective small width hamiltonian dynamics. Note 
also that the introduction of the gluon mass (i$ in the regularization could matter for the lower bound on 
\x — y\ and it could even eliminate the whole contribution when the upper bound of u> meets the lower bound 
of [is- We assume here fig = 0. 

The divergent part in Eq. (3.116) is independent of r and it is easily removed by subtracting 1 from 
exp iqr. Evaluation of the integral leads to the answer that for large r the seagull term produces a logarithmic 
potential of the form 

^ ^Mer) logr , (3.117) 

where a is equal 1 for the radial versor e r along the z-axis and it equals 2 when r is purely transverse. This 
potential is confining. It is also boost invariant. But the rotational asymmetry of the potential raises doubts. 
It suggests that an important piece of physics is missing in the reasoning used to derive it. The obvious 
sources of questions are the mechanism of blocking the effective gluon emissions and absorptions, role of the 
operation R, the role of the nonrelativistic approximation, the size of the quark and gluon masses and the 
strong dependence of the term on the width A. The most urgent question is what happens in higher order 
calculations. 



4. CONCLUSION 

We have defined and illustrated on a few perturbative examples a general method of calculating light-front 
hamiltonians which can be used for the relativistic description of interacting particles. The starting point in 
the calculation is a field theoretic expression for the bare hamiltonian density. This expression is multiply 
divergent in the physically interesting cases. Therefore, the hamiltonian theory requires renormalization. 
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In the renormalization process, one calculates a whole family of effective hamiltonians as functions of the 
width parameter A which determines the range of the effective interactions on the energy scale. 

An effective hamiltonian of a small width A is much different from the initial bare hamiltonian. It 
couples only those states whose masses differ by less than a prescribed amount. Thus, the effective theory 
contains only near-neighbor interactions on the energy scale. No scale is removed in the calculation but 
the correlations between dynamics at significantly different energy scales are integrated out. Therefore, in 
principle, the effective eigenvalue problem can be solved scale by scale using standard techniques for finite 
matrices which describe dynamics at a single scale. 

Our formalism is based on the earlier work on renormalization of hamiltonians from Refs. Q and || 
where the hamiltonians are defined by their matrix elements in a given set of basis states. Wegner has 
developed similar equations for hamiltonian matrix elements in solid state physics. ^] The present approach 
to renormalization of hamiltonians introduces the following features. 

Our similarity transformation is defined in terms of creation and annihilation operators. Consequently, 
calculations of counterterms in perturbation theory can be performed without knowing details of the specific 
Fock states which are needed to evaluate the matrix elements. This is useful because a large number of Fock 
states needs to be considered. The renormalization scheme is free from practical restrictions on the Fock 
space sectors. 

Expressing the effective hamiltonians in terms of the creation and annihilation operators of effective 
particles and showing that the effective interactions are connected is a prerequisite to obtain the cluster 
decomposition property. [ pp| The effective interactions in our approach do not contain disconnected terms. 
The number of creation and annihilation operators in a single term is limited in perturbation theory by 
2 + n{V — 2) where n is the order of a perturbation theory and V is the number of operators in the 
perturbing term. 

The physically motivated assumptions about the model space of effective states included in solving a 
particular problem are introduced after the effective hamiltonian is calculated. The interaction terms in the 
effective hamiltonian contain the similarity factors which diminish the dynamical significance of the Fock 
sectors with numbers of effective particles considerably different from the number of effective particles in the 
dominant sectors. 

The present operator formulation does not introduce spectator dependent interactions, even in the case 
where we include the sums of the invariant masses for incoming and outgoing particles in the similarity 
factors. The sums are useful for estimates of cutoff dependence in perturbation theory. 

The formalism explicitly preserves kinematical symmetries of the light-front frame. The structure of 
counterterms is constrained by these symmetries, including boost invariance. Hence, the number of possible 
terms is greatly limited. Preserving boost invariance is particularly important because it is expected to help 
in understanding the parton model and constituent quark model in QCD, simultaneously. 

It is essential to include the running of the coupling constants in the calculation of the small width 
dynamics. The examples of second order calculations we described in this article do not include the running 
coupling constant effects. Inclusion of these effects requires higher order calculations. 

Wegner 's equation can be adapted to building an operator approach similar to what we described in the 
present article. The initial equation which replaces our Eq. (2.29) when one uses the Wegner generator of 
the similarity transformation is 

= ^ [^-^Wa] • (4-1) 

However, there is little flexibility left in the equation so that the widening of the hamiltonian band is not 
readily available. 

There exists a class of generalized equations for the flow of hamiltonian matrix elements described in 
Ref. JtJ and already studied in a simple numerical model. These equations allow widening of the effective 
hamiltonian matrix at large energies. The generalized equations can also be adapted for the construction of 
the creation and annihilation operator calculus. Namely, 

^ = {F{H 2X }, H X ] . (4.2) 

These equations require detailed definitions of the similarity factors generated by the operation F. 

In summary, the present formalism for renormalization of hamiltonians in the light-front Fock space 
provides a tool for working on a host of theoretical issues in particle dynamics. Second order applications 
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produce boost invariant Yukawa potential, Schrodinger equation for internal bound state dynamics and log- 
arithmically confining quark-anti-quark interaction. However, it remains to be verified if the formalism can 
lead to quantitative improvements in our description of particles. Rotational symmetry and infrared singu- 
larities in gauge theories require further studies. Most urgent are the calculations of effective hamiltonians 
in the third and fourth order perturbation theory. 
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